F\U  COPV  ^ 


AFGL-TR-82-0021 


CJHAVITY  GRADIOMETER  SURVEY  AND  REAL  TIME 
TECHNIQUES  FOR  M PROVING  INERTIAL 
NAVIGATION  SYSTEM  ACCURACY 


!  Daniel  B.  DeBra 
John  V.  Breakwell 
Eugene  Wells 


Stanford  University 
Guidance  &  Control  Laboratory 
Department  of  Aeronautics  and  Astronautics 
Stanford,  California  94305 


Final  Report 

January  1978  -  December  1981 


December  1981 


Approved  for  public  release;  distribution  unlimited 


AIR  FORCE  GEOPHYSICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
HANSCOM  AFB,  MASSACHUSETTS  01731 


82 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (IWifn  Oiic^ Enlereilj 

I  REPORT  DOCUMENTATION  PAGE 


3  Anp  READ  INSTRUCTIONS 

^ _ BEFORE  COMPLETING  FORM 

2.  GOVT  ACCESSION  NO.  3.  RECIPIENT’S  CATALOG  NUMBER 


I.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO.  3.  RECIPIENT’S  CATALOG  NUMBER 

AFGI/-TR-8»-0021  4^-/4//^ 

4.  title  Cand  5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Gravity  Gradlometer  Survey  and  Real  Time  Techniques  Final  Report 
for  Improving  Inertial  Navigation  System  Accuracy  Jan  78  -  Dec  1981 

6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  authorc*; 

Daniel  B.  DeBra 
John  V,  Breakwell 

_ Eugene  Wells _ _ 

9.  PERFORMING  ORGANIZATION  name  ANO  ADDRESS 

Stanford  University 

Guidance  &  Control  Laboratory 

g!g^fo^§T°eg“^|t§0§  Astronautics _ 

It.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Air  Force  Geophysics  Laboratory 
Hanscom  AFB,  MA  31731 

_ Monltot;  James  Hammond/LWG _ 

l4.  monitoring  agency  name  a  ADORESScJ  dlZ/erenl  Itom  ConUatUnt  Otiicc) 


8.  CONTRACT  OR  GRANT  NUMBERr>J 

F19628-78-C-0038 


to.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

63701  B 
3201 03AE 

12.  report  date 

December  1981 

13.  NUMBER  OF  PAGES 

124 

IS.  SECURITY  CLASS,  (ol  this  report) 

Unclassified 


ISa.  DECLASSIFICATION/ downgrading 
SCHEDULE 


1 16.  Distribution  statement  (oi  ihia  Rtpon) 


Approved  for  public  release;  distribution  unlimited. 


t7.  Distribution  statement  (of  tho  mbatroct  ontarod  in  Block  20,  if  different  from  Report) 


I  16.  SUPPLEMENTARY  NOTES 


- iO 


19.  KEY  WORDS  on  r«v«f»o  aldo  it  neceaamry  and  Idonfify  by  block  numbor>  ^ 

gravity  survey,  gravity  gradlometer,  inertial  navigation,  system  accuracy, 
second-order  model,  Heller  model,  cross-track  position  error,  real-time 
models,  covariance  comparisons. 


1 20.  ABSTRACT  (Continue  on  reverae  aide  If  neca«««ry  and  identify  hy  block  number) 


see  back  side 


-  ,  I  « 

1  w  1  w  ;  •  *''4  • 

_ . _ St  ;  '  Ju  a_.a_ 


9D  1  jan*73  1473\i  ^OtTION  of  I  NOV  6S  IS  OBSOLETE 


Unclassified  •  — 

A E CuRlTY  CLASSIFICAT'*"  THIS  PAGE  fWheti  Dmtm  Snfaradt 


Unclassified _ 

security  classification  of  this  PAGEflWiii  Oaf  Enimrad) 


Abstract 


As  inertial  navigation  instrumentation  become!  more  advanced, 

errors  associated  with  imprecise  knowledge  of  the  true  gravity  vector 

\ 

become  significant.  With  the  development  of  the  moving  base  gravity 
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true  gravity  vector  is  available.  This  research  explores  both  real  time 
and  a  priori  schemes  for  using  gradiometer  information  to  reduce  the 
position  and  velocity  error  associated  with  an  Inertial  navigation  system. 
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Chapter  I 


INTRODUCTION 

A.  PROBLEM  STATEMENT 

For  inertial  navigation  a  gyro  stabilized  platform  is  used  to 
spatially  orient  accelerometers  which  measure  specific  force.  Specific 
force  is  the  vector  (difference)  of  the  quantity  to  be  determined,  inertial 
acceleration,  and  gravity.  Therefore  to  Isolate  inertial  acceleration 
the  true  gravity  vector  must  be  found.  It  cannot  be  measured  directly 
from  a  moving  vehicle  because  it  is  inseparable  from  the  effects  of  inertial 
acceleration.  However,  the  gradients  of  the  gravity  field  are  separable 
from  those  of  a  linear  acceleration  field  (zero  if  nonrotating)  thus  making 
it  possible  to  calculate  the  changes  in  gravity.  By  spatially  integrating 
these  changes  on-line,  variations  in  the  true  gravity  vector  may  be  deter¬ 
mined.  If  these  variations  are  used  in  conjunction  with  periodic  platform 
leveling  and  position  fixes,  the  true  gravity  vector  can  be  found.  It  is 
now  a  simple  matter  to  difference  the  true  and  assumed  (from  a  reference 
gravity  model)  gravity  vectors  to  obtain  the  gravity  perturbation  vector. 

The  horizontal  components  of  the  gravity  perturbation  vector  represent  the 
deflections  of  the  vertical  uncertainties  of  which  cause  errors  in  the 
horizontal  inertial  navigation  channels. 

The  problem  addressed  by  this  research  is  to  determine  how  to  use 
the  gradiometer  to  most  effectively  reduce  the  gravity  deflection  error 
contribution  to  the  inertial  navigation  system  (INS)  horizontal  position 
and  velocity  outputs.  This  is  distinct  from  the  greater  question  of  how 
large  the  geodetic  error  Is  when  compoicd  with  the  instrument  errors  in 
the  inertial  navigation  system.  A  discussion  of  this  separate  question 
is  given  by  Hildebrant,  et  al.,  [1]  and  shows,  for  Instance,  that  the 
geodetic  portion  of  the  total  inertial  navigation  velocity  error  budget 
may  be  as  high  as  70%  in  state-of-the-art  systems.  There  are  two  ways 
to  use  the  gradiometer  to  Improve  Inertial  navigation.  Tlie  first  way 
is  to  estimate  in  real  time  the  velocity  and  position  errors  propagating 


in  the  INS  basic  output  loop.  In  reality  this  portion  of  the  problem 
primarily  involves  the  derivation  and  evaluation  of  approximate  low- 
order  inters  for  velocity  and  position  error  estimation.  Approximate 
filters  are  advantageous  because  they  circumvent  the  problems  which 
arise  from  working  with  Inherently  transcendental  Rravity  field  spectra. 

The  second  method  is  to  svirvey  the  area  which  the  INS  carrier  vehicle 
Is  to  traverse.  Using  this  information  the  reference  gravity  model 
contained  in  the  INS  output  loop  is  updated  to  reduce  the  residual  be¬ 
tween  true  and  reference  gravity.  This  involves  updating  the  gravity 
perturbation  power  spectrum  and  subsequently  determining  rms  position  and 
velocity  error  for  survey  parameter  optimization  and  real  time  comparisons. 

Since  navigation  on  the  surface  of  the  earth  is  generally  more 
important  than  vertical  position,  this  research  will  focus  on  only  the 
horizontal  channels  of  a  locally  level  INS,  though  some  of  the  methods 
could  be  applied  to  the  vertical  channel. 


B,  GHADIOMETER  DEVELOPMENT  AND  PREVIOUS  WORK 

Moving  base  gravity  gradiometers  have  been  under  development  for  the 
past  13  years.  This  development,  done  under  the  auspices  of  the  Air  Force, 
Navy  and  the  Defense  Mapping  Agency,  resulted  in  two  workable  gradiome— 
ters;  one  by  the  Bell  Laboratories  [2],  and  one  by  the  Draper  Labs  [3].  The 
Bell  instrument  was  chosen  by  the  Navy  as  the  most  advantageous 
design  for  their  application,  and  will  be  used  when  reference 
to  hardware  is  necessary.  The  Bell  instrument  is  commonly  called  a 
"rotating  gravity  gradiometer"  (RGG)  because  the  measurement  unit  is 
actually  rotated  to  separate  spectrally  the  gradient  signals  from  out¬ 
side  disturbances. 

The  gradiometer  research  done  at  Stanford  began  in  1972 
under  the  auspices  of  the  U.S.  Air  Force  Geophysics  Laboratory.  Pelka's 
dissertation  f^"]  dealt  with  improving  the  accuracy  of  the  Hughes' 

(prototype)  gradiometer  using  on-line  parameter  estimation  and  control 
techniques.  Subsequent  studies  by  DeBra,  Breakwell,  Schaechter,  and 


Ktirosnki  deiilL  with  thr  clfvelopiiicMiL  of  K*'avity  iiiotlfls  used  to  <lc(cr- 
mine  the  performance  of  the  fcradiometer  in  mappins  the  earth's  nt'avity 
field  [5],  The  culmination  of  this  work  describes  the  use  of  a  gravity 
gradlometer  aboard  a  horizontally  moving  vehicle  to  estimate  gravity 
deflection  and  anomaly  [6].  Gradlometer  derived  improvements  of  inertial 
navigation  and  gravity  survey  work  are  presented  in  this  dissertation. 


C,  NEW  RESULTS 

Whereas  most  previous  research  on  gradiometers  has  dealt  witli  the 
instruments  themselves,  the  research  reported  here  deals  with  the 
processing  of  gradlometer  data  both  in  real  time  and  after  the  fact 
(the  latter  being  called  the  "survey"  use).  No  other  investigators 
have  addressed  these  aspects  of  the  problem  as  they  have  been  addressed 
here. 


Specific  contributions  are: 

1.  The  development  of  a  method  for  deriving  approximate  Wiener 
filters  for  estimation  of  INS  position  and  velocity  error, 

2.  The  statistical  evaluation  of  these  filters  to  determine  rela¬ 
tive  and  absolute  optimality; 

3.  The  statistical  evaluation  of  various  survey  schemes  for  update 
of  the  reference  gravity  model; 

4  Development  of  an  improved  version  of  Heller’s  gravity  model 

to  portray  more  accurately  the  perturbation  potential  in  light 
of  the  latest  gravity  data  available. 
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1).  THESIS  OUTLINE 

The  full  inertial  navigation  system  error  equations  and  the  assump¬ 
tions  used  to  simplify  them  for  use  in  this  researcli  are  presented  in 
Chapter  II. 

Chapter  III  presents  three  statistical  p:ravlty  perturbation  models 
used  in  the  covariance  analyses  presented  in  Chapters  V  and  VII,  The 
first  model  is  rational  (the  reason  for  its  use)  but  somewhat  inaccurate. 
The  second  and  third  models  are  transcendental  but  are  representative 
of  the  most  accurate  gravity  perturbation  models  available. 

In  Chapter  IV  the  basic  theory  used  in  deriving  approximate  real 
time  filters  is  shown,  along  with  examples,  for  a  wide  spectrum  of  meas¬ 
urement  combinations. 

Several  filters  obtained  in  Chapter  IV  are  evaluated  and  compared 
by  covariance  methods  in  Chapter  V. 

In  Chapter  VI  the  theory  used  in  analyzing  the  effect  of  a  gravity 
survey  on  INS  performance  is  developed. 

In  Chapter  VII  covariance  data  for  the  survey  scheme  is  computed 
and  compared  with  data  obtained  from  the  real  time  application.  The 
ultimate  goal  of  this  comparison  Is  to  determine  which  gradlometer 
application  is  most  desirable  under  a  given  set  of  conditions. 

In  Chapter  VIII  the  conclusions  and  future  recommendations  of  this 
research  are  presented. 


Chapter  II 

INERTIAL  NAVIGATION  ERROR  EQUATIONS 


Before  looking  at  the  error  equations  associated  with  conventional 
Inertial  navigation  systems  (INSs)  it  is  instructive  to  examine  the 
basic  navigation  loop  used  throughout  this  research.  Figure  II-l  is 
a  schematic  of  this  loop.  Note  that  some  type  of  external  velocity 
measurement  (v"*)  is  assumed  for  damping  ol'  Schuler  oscillations. 

A.  FULL  INERTIAL  NAVIGATION  SYSTEM  ERROR 
EQUATICTJS 

The  error  equations  relating  the  gravity  perturbations  to  INS 
position  and  velocity  error  have  been  developed  in  detail  by  Heller  [7]. 
These  equations  assume  an  INS  mechanized  in  a  local— level  wander— azimuth 
frame  and  coordlnatized  in  a  local— level  north-pointing  frame.  The 
full  position  and  velocity  error  equations  for  the  horizontal  channels 
of  an  INS  are  given  below  [from  Ref.  7]: 

%  =  »*N  ■  ^’^E  Vz  -  '  «N 

V-; 

( 2il  sin  L  +  —  tan  L)  6v 

*  - "%  -  "V.  *  "Ve  * 

V  V 

+  —  tan  L)f.V  +  (2W  cos  L  +  — )i'V 
Ry  N  H,.  z 

-  -  kV^t^  +  kVj^tlr^  +  k:>v"  ; 
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FIG,  II-l  SCHEMATIC  OF  BASIC  NAVIGATION  LOOP 
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JB.  SIMPLIFYING  ASSUMPTIONS  AND  RESULTING 
ERROR  EQUATIONS 

The  error  equations  are  linear  thus  aHowins  different  error  con¬ 
tributors  to  be  superposed.  In  order  to  focus  on  gravity  induced  errors 
the  nongravltational  error  sources  in  (2.1)  to  (2.4)  (i.e.,  accelerometer 
bias,  p,,  platform  misalignment,  etc.)  may  l)e  set  aside  [iT.  Further¬ 

more,  the  north  and  east  equations  may  be  decoupled  without  limitation  to  re¬ 
search  objectives.  Therefore  the  following  assumptions  have  been  made. 

1)  The  carrier  vehicle  travels  due  north*; 

2)  There  are  no  platform  misalignment  or  accelerometer  bias  errors; 

3)  The  carrier  vehicle  altitude  is  fixed; 

4)  The  earth  Is  nonrotating. 

With  these  assumptions  the  horizontal  channel  error  equations  become 
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This  assumption  is  simply  a  manifestation  of  the  particular  north-east 
down  coordinate  system  chosen  in  (2.1)  through  (2.4).  Equations  (2.5) 
through  (2,8)  would  remain  the  same  had  another  local-level  coordinate 
system  been  selected  to  coincide  with  the  In-track  (x)  and  cross-track 
(y)  directions  of  a  vehicle  traveling  other  than  north.  Thus  the 
analysis  done  here  applies  to  a  carrier  vehicle  traveling  In  any  direc¬ 
tion  over  a  non-rotating  earth. 
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Note  that  since  the  carrier  vehicle  travels  due  north  this  direc¬ 
tion  will  henceforth  be  denoted  irv-track  (x  subscript)  and  the  easterly 
direction  denoted  cross— track  (y  subscript). 


Chapter  III 


STATISTICAL  C.KAVITY  PEIITURBATION  MODELS 

A  stochastic  study  of  INS  position  and  velocity  error  requires  that 
the  gravity  perturbation  terms  in  Eqs.  (2.5)  and  (2,6)  be  statistically 
modeled.  Three  models  were  used  in  this  study,  Tlie  first  is  a  relative! 
inaccurate  rational  model  used  solely  for  Kalman  filter  comparisons,  as 
a  rational  model  is  necessary  for  a  Kalman  realization.  The  last  two 
models  have  the  same  basic  form,  the  only  difference  being  a  slight  para¬ 
meter  modification  in  the  final  version  to  fit  observed  data  more  closely 
These  models  are  transcendental  but  portray  more  accurately  the  earth's 
actual  gravity  field.  All  throe  of  these  mo<)els  were  developed  by  the 
processing  of  empirical  data  such  as  that  collected  by  Rice  [8]. 

_A.  SECOND-ORDER  MODEL 

The  details  of  this  model  are  fully  detailed  by  Kasper  [9]  and  reit¬ 
erated  by  Heller  [10].  To  summarize  briefly,  it  is  a  second-order  model 
consisting  of  zero  mean  white  noise  put  through  a  shaping  filter  as  shown 
in  Fig.  III-l. 


FIG.  III-l  SECOND-ORDER  GRAVITY  PERTURBATION  MODEL 

The  Markov  parameter  P  is  defined  as 

p  =  V^/d*  .  (3.1) 

where  d  is  the  vertical  deflection  ensemble  chnractc'ristlc  distance. 


The  constant  d  in  Kasper's  model  Is  related  to  the  more  commonly 
known  correlation  dlstanc<*  (D)  by  the  relation:  D  =  2.l‘16d. 
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This  model  may  be  iise<l  in  either  the  in— track  or  cross-track  direc¬ 
tions,  However,  the  Kalman  filter  coi'iparisons  need  only  be  done  cross- 
track  to  demonstrate  the  accuracy  of  the  approximate  filters.  The  spec¬ 
tral  density  of  the  white  noise  ,  Is  given  as 
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A  state  space  description  of  this  model,  required  in  Ch,  V,  may  be  written 
as 
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where  is  a  Markov  Process  state  variable. 


BASELINE  HELLER  MODEL 

This  statistical  model  was  developed  by  Heller  [11]  and  descrilies 
perturbations  in  the  earth's  surface  potential  away  from  an  ellipsoidal 
gcoid  as  fictitious  white  noise  potential  fluctuations  restricted  to  a 
few  thin  layers  below  the  earth’s  surface.  The  model,  as  shown  below, 
is  defined  by  tlic  power  spectral  density  (PSD)  of  surface  potential 
fluctuations 
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The  values  for  the  depths  (U^*s)  and  PSDs  of  the  three 

layers  of  interest  given  by  Heller  are  listed  in  Table  III-l. 


Table  III-l 

BASELINE  HELLER  MODEL  PA[{AMETERS 


To  determine  the  relationship  between  surface  potential  fluctuations 
and  potential  fluctuations  at  an  altitude  z  above  the  earth’s  surface 
it  facilitates  analysis  to  work  with  a  "flat  ear'^h"  approximation  [12] 
to  Poisson's  upward  continuation  formula  for  the  earth's  potential, 
Poisson's  upward  continuation  formula  is  given  by 
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(3.5) 


where  the  direction  of  the  vector  r  is  determined  by  the  (r,  ib)  spher¬ 
ical  coordinate  set,  and  the  direction  of  is  determined  by  the  set 

(R^,  fl',  (b').  The  integration  of  (3.5)  is  done  over  the  entire  earth's 
surface.  The  details  of  the  derivation  of  (3.5)  are  given  in  App.  A. 

Making  the  "flat  earth"  approximation  (z  <<  R^^ )  to  (3.5) 
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where  (E)  denotes  convolution  in  two  <liniensions .  Taking  the  fourier 
transform  of  yields 
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Tills  allows  (3.4)  to  ho  rewritten  as 
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Uecause  Heller  has  modeled  the  potential  perturbations  it  is  a 
simple  matter  to  obtain  from  (3,8)  the  two  dimensional  PSD  expressions 
for  any  of  the  gravity  or  gravity  gradient  perturbation  components.  In 
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where  all  w's  Indicate  spatial  frequency  variation. 


The  spectral-density  matrix  along  the  line  y  -  0  is  required 
because  the  carrier  vehicle  t .s ,  by  definition,  always  traveling  in-track. 
This  matrix  is  obtained  by  performing  the  following  operation: 
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(3.9) 
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The  results  of  this  operation  are  again  shown  by  Breakwell  and  Kurosaki  [13] 
as 
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where 

Kj(ni)  =  modified  Bessel  function  of  first  order 
=  normalized  frequency  variable  (0^  = 

'  =  differentiation  with  respect  to  the  variable  0^^  . 

g 

Although  in  [m]  the  designated  states  are  g  ,  g  ,  and  g  , 
i*l'  X  y  z 

the  PSDs  of  the  In-track  spatial  derivatives  of  these  quantities  V 
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r  may  be  found  through  the  application  of  Laplace's  equation 
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The  zero  correlations  in  (3.10)  are  Intere.stinK  because  they  show 

wliich  (iradiont  pertnrl>ation  components  iiiav  be  used  to  estimate  and  a 

X  y 

in  Eqs .  (2.5)  and  (2.0).  The  zero  correlation  between  s  ,  P  and 

y  X  y 

the  other  states  indicates  that  the  optimal  estimate  of  should 

involve  only  the  measurements  T  and  r  .  The  optimal  estimate  of 

xy  yz 

K  ,  on  the  other  haml ,  should  use  only  T  ,  T  ,  T  (or  T  ). 

X  XX  xz  zz  yy 

Spectral  separation  considerations  to  be  explained  in  Ch.  V  dictate 
that  a  low  frequency  approximation  to  (3,10)  be  available.  Such  an  approxi¬ 
mation  was  developed  by  Breakwell  and  Kurosaki  [13]  and  is  shown  below. 
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IMPItOVEI)  HELLER  MODEL 

In  the  course  of  this  research  it  became  necessary  to  improve 
Heller’s  baselltie  model  with  the  latest  pfravlty  data  available  in  an 
attempt  to  shorten  the  characteristic  times  of  certain  low  frequency 
poles,  to  l)e  introduced  in  the  next  cliapter.  The  process  and  result 
of  this  improvement  will  be  di.scussed  here. 
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The  data  used  to  improve  the  baseline  Heller  model  (BHM)  was  the 
Goddard  Earth  Model  (GEM)  10  data  published  by  NASA  [i'll.  In  the  GEM  10 
studies  done  by  the  NASA  the  earth’s  potential  was  described  in  the  stan¬ 
dard  spherical  harmonic  form  Kiven  below; 
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where  the  C«  's  and  Sd  's  determined  from  satellite  data  actually  define 
•*  m  *  m 

GEM  10. 


If  it  is  assumed  that  the  gravity  model  used  in  an  INS  is  at  least 
as  accurate  as  GEM  10  then  the  rms  potential  coefficient  error  estimates 
for  GEM  10  taken  from  Ref.  14  (p.  05)  should  equal  the  rms  potential  coef¬ 
ficient  variation  (a-  )  of  the  Heller  model.  The  relation  defining  the 
mean-squared  potential  coefficient  variation  for  both  the  spherical  and 
Heller  models  is  derived  in  App.  B  and  summarized  as 
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By  substituting  the  GEM  10  rms  coefficient  error  estimates  in  for  C 

I'm 

and  S^^  in  (3.15)  the  PSD  of  the  deepest  layer  of  Heller's  model  may 
be  adjusted  so  that  (3.15)  is  roughly  satisfied  at  low  frequency. 

Recalling  that  the  reason  for  this  exercise  is  to  adjust  low  frequency 
pole  characteristic  times,  only  need  be  adjusted  as  this  value  has 

the  largest  effect  on  low  frequency  behavior  while  its  influence  dies  out 
at  high  frequency.  A  satisfactory  match  was  obtained  by  reducing 
by  a  factor  of  10.  Figure  III-2  .shows  the  results  of  this  change.  Thus 
the  IHM  (improved  Heller  model)  parameters  are  the  same  as  tliose  given 
in  Table  III-l  except  that 

=  1.16  X  10^  km^/sec"*  .  (3.16) 
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The  improved  Heller  model  ( IHM)  was  a  necessary  development  In  the 
course  of  this  research  as  explained  above.  To  alleviate  any  confusion 
in  deciding  which  Heller  model  was  used  in  the  comparisons  to  come, 
great  care  has  been  taken  to  clearly  mark  each  case. 

Thus  far  in  the  discussion  all  of  the  power  spectral  densities  have 
been  based  on  spatial  as  opposed  to  time  variation.  In  the  following 
chapters,  PSDs  based  on  time  will  bo  required.  The  conversion  from 
spatial  to  time  variation  for  a  carrier  vehicle  with  constant  velocity 
simply  involves  the  carrier  speed  as  follows 

r  X  V  .  (3.17) 

Note  also  the  spatial— to-tiine  frequency  variable  transformation 

t  s 

w  ^  ui  X  V  (3.18) 

which  will  be  used  for  consistency  in  writing  various  PSD  expressions. 


Chapter  IV 


APPROXIMATE  REAL  TIME  FILTERS 


In  this  chapter  the  basic  theory  for  deriving  filters  for  the  estima¬ 
tion  of  INS  position  and  velocity  error  is  presented.  Gradiometer  outputs, 
and  in  some  cases  the  externally  obtained  vehicle  velocity  measurement, 
are  used  as  inputs  to  these  filters.  Figure  IV-1  shows  a  schematic  of  the 
real  time  gradiometer  application  scheme  which  is  simply  added  on  to  the 
basic  navigation  loop  shown  in  Fig.  II— 1, 

The  scope  of  this  thesis  permitted  at  most  two  gradiometer  measure¬ 
ments  In-track  or  cross-track.  The  output  signal  of  a  single  Bell  RGG 
with  its  spin  axis  vertical  is 

-r(r  -  r  )  sin  2flt  +  y  cos  2nt^ 

2  '  yy  '  XX  'xy 

where  y  is  the  gravity  gradient  tensor  prior  to  differencing  with  the  ref¬ 
erence  gravity  gradient,  Q.  is  the  gradiometer  rotation  rate,  and  t  is  time. 

After  demodulation  to  separate  the  sine  and  cosine  terms  and  differencing 
the  gradient  components  with  the  reference  values,  the  following  outputs 
are  obtained 
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where  r  is  the  difference  gradient. 


The  navigation  axes  x  and  y  to  which  the  sensor  is  aligned  are  not 

necessarily  principal  axes  of  the  gravity  gradient.  In  the  plane  of 

rotation  of  the  sensor  the  principal  axes  are  at  an  angle  with  respect 

to  the  X  and  y  axes.  If  a  pha.se  angle  is  introduced  in  the  equation 

indicated  by  this  footnote  it  would  indicate  the  orientation  of  the 

principal  axes  when  the  cross-gradient  component  is  zero.  Tiuis  the 

gradient  can  also  be  represented  by  an  angle  and  the  maximum  difference 

(y  -  y  )  in  the  plane  of  rotation, 

'  yy  XX 

t  The  outputs  obtainable  from  a  single  spliericnl  gradiometer  manufactured 
by  the  Draper  I.nb.  are  *lifferent  in  form  from  those  of  the  Bell  RGG. 

It  is  an  easy  matter,  however,  t€>  alter  I  he  anatysis  given  in  tills 

thesis  to  accommodate  tlie  Draper  gradiometer.  An  overview  of  the 

operation  and  outputs  of  the  spherical  gradiometer  is  given  by  Trageser f 15 1 . 
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FIG,  IV-1  SCHEMATIC  OF  REAL  tIME  FILTER  APPLICATION 
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A  cyclic  permutation  of  these  subscripts  {jives  the  outputs  for  instruments 

with  spin  axes  aligned  in-track  and  cross-track  respectively.  Ilecall 

from  Ch.  Ill  that  only  the  cross-gradients  r  and  r  can  be  used  as 

xy  yz 

measurements  cross-track  and  only  the  iit-line  gradients  (r  ,  f  P  ) 

XX  yy  7.7, 

and  can  be  used  as  measurements  in-track.  Thus  to  obtain  two  gradi- 

ometer  measurements  cross— track,  two  instruments  are  required  with  spin 

axes  in-track  and  vertical.  The  in-line  gradient  combinations  from  these 

same  two  instruments  are  used  for  in-track  estimation.  This  allows  for  a 

full  evaluation  of  how  two  gradtometers  can  improve  horizontal  navigation. 

The  best  choice  for  a  single  cross— track  measurement  is  P  because  it 

xy 

is  more  highly  correlated  with  g  than  is  P  ,  For  consistency  the  in- 

y 

line  output  of  this  same  gradiometer  is  used  as  a  single  in-track  measure¬ 


ment  . 


Table  IV-1  lists  each  of  the  measurement  combinations  used. 


Table  IV-1 

CROSS-TRACK  AND  IN-TRACK  MEASUREMENT  COMBINATIONS 


Measurement 

Combination 

Code 


Pi r  fcrcncc 
Gradient 
Con«poncnl.s 


pr  *  Ir  1 

'■XX  /  r. 


r  ♦Jr 

XX  *  /7 


4  r 

•  XX  r.r. 
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The  significance  of  each  is  given  helow. 


CT-1,  This  combination  corresponds  to  using  the  cross-gradient  out¬ 
put  of  one  Bell  rotating  gravity  gradiometer  (RQG)  witd  its  spin 
axis  mounted  vertically.  Because  this  tensor  cppiponent  is  the 
in-track  spatial  derivative  of  the  cross-track  deflection 
it  is  the  tensor  component  most  highly  correlated  with  g 

y’ 

CT-2.  To  obtain  this  set  of  measurements  two  ItGGs  are  required. 

the  spin  axes  of  these  instruments  point  yerticsily  and  In- 
track  respectively.  The  improvement  over  the  CT-1  combina¬ 
tion  is  small, 

IT-1.  Because  tlio  in-line  gradients  come  from  the  Bell  RGG  linearly 

combined  (e.g.,  for  spin  axis  vertical  i(r  -  T  ))  V 

XX  yy  ’  XX  ’ 

the  most  highly  correlated  component  wjth  g  is  npt  directly 
obtainable.  However,  for  comparison  purposes  this  fictitious 
case  of  a  direct  measurement  is  included. 

IT-2,  Tills  case  illustrates  the  situation  where  the  direct  in-line 
composite  oqtput  (modified  by  Laplace's  equation)  of  one 
vertical  RGG  Is  used  as  a  measqrement .  This  is  inferior 
to  the  IT-1  case. 

IT-3.  In  this  situation  the  direct  composite  outputs  of  two  RGGs 

arc  used  where  the  instrument  spin  axes  are  vertical  and  in- 
tracR.  This  combination  is  comparable  in  information  to  the 
IT-1  case. 

_A.  VELOCITY  ERROR  ESTIMATION  THEORY 

1  Excluding  the  External  Vclo<'lty  Measurement 

To  obtain  the  measurement  combinations  .just  discussed  the  basic 
difference  grailicnt  vector  required  is 


t 


xy 


yz 


The  measurements  may  then  be  written  as 


(4.1) 


myl  mx5  5x1  mxl 

z  ^  II  f  +  i 


(4.2) 


where  e  Is  the  white  measurement  noise  vector. 

Let  the  estimated  velocity  error  lx? 

IXm  )i^l 

Sv  (s)  =  ijf(s)  z  =  ilf(s)[Hr  I-  g]  . 


(4.3) 


where  'li(s)  is  the  velocity  error  filter  expre.ssed  in  the  Laplace  domain. 


This  vector  requires  only  these  five  components  because  any  sinule 
channel  output  of  a  Bell  RCKJ  with  axes  oriented  alons  the  x,  y,  z 
directions  can  be  reproduced  with  these  elements  and  the  proper  H 
matrix.  If  all  the  elements  of  the  gravity  perturbation  tensor  arc 
required  three  orthoRonal  IlGGs  are  required.  These  instruments  must 
be  mounted  so  that  a  slKiii ficant  component  of  gravity  is  perpendicular 
to  the  spin  axis  of  each  instrument.  This  is  necessary  for  proper 
operation  of  the  scale  factor  calibration  loops. 


Note  here  that  (4.3)  is  written  with  subscript  x  Indicating  the  in¬ 
track  (iirection.  nocause  the  basic  INS  error  equations  have  the  same 
iorm  IT  or  CT  the  velocity  error  estimation  equations  will  be  written  IT 
but  are  immediately  applicable  CT  simply  by  chanKlng  the  x  subscripts 
to  y.  Taking;  the  Laplace  transform  of  (2.5)  and  (2.7)  and  solving  for 
fSV  gives 

X 

m 

sg  skRV 

fSV  (s)  =  - -  + - -  (4.4) 

(+)  (+) 

who  re 

(+)  -  s^  +  ks  ^  (4.5) 

s 

is  called  the  positive  Schuler  factor  and  conversely,  (-)  will  be  used 

2  2 

to  denote  the  negative  Schuler  factor  (s  -  ks  +  .  Subtracting  (4,4) 

from  (4.3)  the  estimate  error  is 


sg  sk6v"' 

5v^(s)  =  »|f(s)Hr  +  ilf(s)e  -  - 


(4.6) 


The  PSD  of  av  (s)  is  then 

X 


(b)s 

g 

s)  +  t(-s)R\lr'^(s)  - 

X  n 


(4.7) 


2  ‘^^^2 


(  s)  (s)hV(s) 


VT 


TT 


mxm 

where  II  is  the  PSD  matrix  of  the  white  noise  measurement  error  vector 
(e).  The  object  of  this  exercise  is  to  minimize  t*'-'  mean-squared  estima¬ 
tion  error  which  may  be  expressed  as 


X 


ds  . 
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1 

2n.i 


(4.8) 


^ ^  ®  miniinum  its  variation  due  to  any  chantje  in  \lr(s) 


5V^ 

should  be  positive.  SubstitutinK  (4.7)  into  (4.8)  and  determining  the 
2 

change  in  cr~y  due  to  a  change  in  ,ir(.s)  gives 

jco  (  H':fg^(s)s 


(s)h’’  +  r1 

.  rr 


s(t  _,j,(s)ir 


(-) 


(s) 


A>1;(-s)|^H<1>I^(s)h^  +  rJ  ii>lf'^(s)  ds 


(4.9) 


Using  contour  Instead  oi  line  integrals  to  evaluate  the  first  two 
terms  in  (4.9)  yields 


.  2 
Aa..- 
oV 


s) 


Hil,  (s)s 
^  '’x 


T!T 


ds 


rS 


^  \  jt{-s)[H<J>^^.j,(s)H’’+R] 


s'lg  i’T(s)n'' 


H - -  |ai1;'^(s)  ds 


+  ^  f  ^'I'(-s)[h'L_.j.(s)h'^  +  r]a«|,'^(s)  ds 

-iOC  ^  ^ 


(4.10) 


where  the  contours  and  are  illustrated  in  Fig.  IV-2  and  enclose 

the  left  and  right  half  complex  plane  respectively.  With  these  con¬ 
tours  it  is  easy  to  see,  using  the  Cauchy  residue  theorem,  that  if 


1) 


s)H^  I  uj  t''^(s)  - 


(s)s 


(+) 


(4.11) 


has  no  left-lialf  plane  poles; 


2) 


The  integrands  In  (4.10) 


converge  properly. 


(4.12) 


Then  the  contributions  of  the  contour  integrals  in  (4,10)  is  zero, 

2 

This  means  that  ACT^  is  positive  since  the  last  integral  in  (4.10) 
is  always  positive  ana  the  minimization  of  o‘~  is  assured. 

Note  from  (4,11)  that  the  poles  of  \|f(s)  iruist  be  determined  either 
by  setting  (+)  -  0  or  by  setting 

*^1  =  0- 

rr 

Because  of  the  inherent  sensitivity  of  the  INS  to  Schuler  frequency 
disturbances  it  will  be  assumed  that  the  dominant  poles  come  from  the 
Schuler  factor  (+).  In  addition,  convergence  in  (4.12)  requires  that  the 
numerator  of  ilr(s)  be  a  well  selected  constant(s).  Therefore  assume 

Ixm 

IXm  - 

t(s)  =  .  (4.13) 

As  promised  \lf(s)  is  of  relatively  low  order.  The  approximate  nature  of 
if\s)  comes  from  the  "Schuler  dominance"  assumption  which  is  made  to  de¬ 
termine  a.  The  idea  is  that  due  to  INS  Schuler  sensitivity  only  the 

Schuler  frequency  portion  of  the  estimate  error  PSb  is  significant  in 
2 

contributing  to  .  This  is  to  say  that  the  estimate  error  PSD  is 

fiVjj 

highly  peaked  at  or  near  Schuler  frequency.  With  this  assumption  condition 

~  t 

(4,11)  is  enforced  only  at  s  to  give 


a’’  =  +  R 

One  way  to  indirectly  show  the  oplimnlily  ol'  I  lie  filler  .jiisl  de¬ 
rived  is  to  prove  the  accuracy  of  the  Schuler  dominance  assumption. 

This  may  be  done  by  substituting  (4.13)  and  (4.7)  into  (4.H)  where  the 

NI 

quadrature  is  performed  numerically  to  give  the  rms  estimate  error  rr~  . 

X 

This  value  can  then  be  compareil  lo  the  analytically  derived  rms  estimate 
AF 

error  obtained  by  assuining  Schuler  ilomlnance.  The  degree  to 

which  agrees  with  rr~  is  a  flirect  indication  of  the  validity 

^  X  X 


(4.14) 
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ol'  assuiiiiiit'  Schuler  dominance  and  an  indirect  indicator  of  filter  optimality. 
AK 

To  dej'ive  Iwiiln  by  writing  only  those  parts  of  (4.7)  which  will 

''x 

liave  left-hall'  plane  poles  if  coiulLtlon  (4.11)  is  assumed  satisfied  across 
the  frequency  spectrum.  Substituti  nj;  (4.13)  into  he  estimate  error  PSD 
thus  obtained  makes  it  possible  to  write  the  mean  squared  estimation 
erroi'  as 


1  /‘F 

sTi  t  IfFir 


2,  2 

s  k  q., 


.(s)hV 


(-*-)(-)  (  +  )  (-) 


^  ds  .  (4.1 


The  IntejfJ'al  in  (4.15)  may  easily  be  evaluated  usinix  the  Cauchy  residue 
theorem  if  the  Schuler  portion  of  the  estimate  error  spectrum  is  again 
assumed  to  be  dominant.  Under  this  premise  and  substituting  (4.14) 
into  (4.15)  tlic  mean  squared  estimate  error  becomes 


C  2  2  2 

•  1  Bw  +  u)  k  q_ 
s  s  2 


(+)  (-) 


(4.16) 


where  B  is  tiie  a  posteriori  gravity  error  spectrum  defined  as 


(  j(,)  )  -  !■  _ 

g  1 


(4.17) 


Then  the  f'au'  liy  residue  evaluation  of  (4.16)  yields 


2  I.  '“•2 

■  2g^—  • 

X 


(4.18) 


Includirif^  the  External  Vclofity  Measuremeiil 


Krom  Section  A-1  recall  that  tlie  coeliicient  of  [\lr(-s)j  il)(-s)]  in 

2 

(•1.2S)  must  have  no  left-half  plane  poles  if  is  to  be  a  minimum. 

Kor  this  to  1k‘  true  the  clenominators  of  i|f(s)  and  'h(s)  must  be  the  left- 
half  plane  roots  which  result  from  settinH:  the  determinant  of  [W^] 
equal  to  x.ero.  Anticipate  that  the  product  of  these  roots  results  in  a 
second  <lenrco  polynomial  of  the  form 

2  2 

(+')  -  s  +  k’s  f  Wg  (4.25) 

where  k’  is  a  sliifted  value  of  k  used  in  the  filter.  To  determine 
the  value  of  k'  tlie  (letenninant  of  (4.24)  is  set  equal  to  zero 
t  he  renion  wlu'i'e 


2 

S  4 


k’s 


2 


+  w 

s 


0 


(4.26) 


or 


2  2 
>s  ■»  lt> 
s 


-k's 


(4.27) 


for  the  determinant  of  [Wj]  to  equal  zero  in  the  Schuler  frequency  region 
the  following  condition  must  hold 


-1 


,1-  /  ■  w.T  ^  t 

X  )H  +  Wg  .  ( jw  ) 

‘  ‘’X  ® 

where  (4.28)  is  obtained  from  reference  to  the  identity 


(4.28) 


Q.. 


lqI  % 


q;‘ 


Lo 


.-I 


(4.29) 
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where  ,  Q^,  represent  the  corresponding  quantities  in  ["^ji¬ 

lt  is  clear  that  if  the  determinant  of  (4.29)  is  equal  to  'oro  tlien 

Substituting  (4.27)  into  (4.28)  yields  the  followinj;  value  I'or  k' 


(4.31) 


Continuinjj  with  the  filter  derivations,  Cor  the  coefficient  ol' 
[i|f(-s);  <t>(s)]  in  (4.23)  to  have  no  left-half  plane  poles,  the  (+) 
factors  in  the  second  column  of  must  be  cancelled  by  the  numerator 

of  ❖(s).  Therefore  assume 


4-  =  - 


(4.32) 


Substituting  (4,32)  into  (4.21)  gives 


sg  sk*oV 

5V  =  t(s)Hf  +  ilf(s)c - - - - 

(  +  ’)  (+') 


(4.33) 


which  is  identical  to  (4.6)  except  that  k'  is  substituted  for  k.  Krom 
this  point  on  the  analysis  is  the  same  as  in  Sect.  A— 1  resulting  in  the 
following  equations 


\|f  (s)  = 


=  +  r]  (.jOg) 


1X1 


k’  -  k)s 


(4.34) 


(4.35) 


(4.36) 


The  analytic  and  numerically  integrated  estimate  error  expressions 
are  also  the  same  as  those  in  Sect,  A-1  except  for  the  change  in  k.  The 
analytic  expression  may  be  written 


(  Al'  ^2 

X 


Note  that  this  is  the  same 
urement  provided  that  k 
error . 


„  k»  q  _ 

u  _ ^  n — I 

~  2k*  2  ~  V^‘*2 


(4.37) 


r  AF  ,2 

(<r~  )  without  the  external  velocity  meas- 

is  optimized  for  the  mean  squared  estimation 


B.  POSITION  ERROR  ESTIMATION  THE(»Y 

In  contrast  to  velocity  error  estimation,  position  error  estimation 
requires  that  tlie  low  frequency,  as  well  as  the  Schuler  frequency,  por¬ 
tion  of  the  estimate  error  spectrum  l)e  considered.  This  is  a  consequence 
of  two  facts.  First  the  high  frequency  content  of  the  spatial  gravity 
perturbation  spectrum  falls  off  quite  rapidly.  T'lis  is  described  in  the 
Heller  model  as  an  exponential  decay.  Tiiis  results  in  the  error  contri¬ 
bution  from  the  narrow  bandwidth  centered  at  Schuler  frequency  becoinlng 
ever  smaller,  for  <lecreasing  carrier  vehicle  velocity,  when  compared  to 
the  low  frequency  spectral  content.  The  second  significant  observation 
is  that  the  gravity  perturbation  spectrum  .just  described  is  integrated 
only  once  to  give  velocity  error  but  it  is  integrated  twice  in  the  posi¬ 
tion  case.  The  low  frequency  error  contribution  is  therefore  much  more 
important  to  position  error  estimation.  This  fact  is  reflected  in  the 
following  assumed  form  for  the  gradiometer  measurement  position  error 
filters . 

ilr(s)  =  \lf,,(s)  +ilf/s)  (4.38) 

where  the  s  and  I  indicate  Schuler  and  low  frequency  portions  of  the 
filter.  The  transfer  functions  i^g(s)  and  '!f£(s)  may  always  be  deter¬ 
mined  1  lulepeiKient  1  y  of  one  another  if  certain  logical  assumptions  are 
made.  These  assumptions,  to  lx?  discussed  in  the  next  section,  concern 
the  relative  size  at  low  frequency  of  terms  contributing  to  the  estimate 
error. 
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B-1  Schuler  Frequency  Filters  (In-track  or  Cross-Track) 

As  with  velocity  error  estimation,  the  Schuler  portion  ol  the  posi¬ 
tion  estimation  equations  will  bo  written  in-track  but  are  iinmefliatel  y 
applicable  cross-track  by  chan^lnK  the  x  subscripts  to  y.  To  estimate 
position  error,  first  without  the  external  velocity  measurement,  hefjin 
as  in  the  velocity  case  with 

^  -  nr  4  K  (4.39) 

and 

W{^(s)  :  ,If(s)z  .  (4,40) 


The  estimate  error  and  estimate  error  power  spectral  density  can  be 
written  as 


m 


?SR  (s)  =  ^lf(.s)H^  4  ilf(s)fi  -  - - 

X  (4) 


(4.41) 


and 


Alf(-s)  _(s)H^  4  Kli|f^(s)  +  - ~7~ 

pp‘  (4)(-) 


(4.42) 


\|f(-s)nil>^  (s)  (J)*"  _._(s)H'^,y^(s) 


(4)(-) 


(4) 


(-) 


FollowlnK  the  same  reasonlnp;  as  in  velocity  error  estimation,  to 
2 


minimize 


Il'T^  (f^> 
^'x 

(s)!!*^  4  H'',/(s)  -  - 

TT  (4) 


(4.43) 
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i 

I 


must  have  no  left— half  plane  poles.  Kocusint;  attention  now  on  Schuler 
frequency  by  assuming  that  i*’  very  small  compared  to  (.jw  ) 

(■1.1S)  may  1h'  fowi'itt<'n  as 

Il4  (s) 

Tg 

t  X  X 

_ ^(s)H  +  ll'],i,^(s) - .  (4,44) 

IT  (+) 


The  position  error  estimation  equivalent  of  (4.12)  dictates  that  the 
numerator  of  %^(s)  be  at  most  linear.  Therefore  a  first-order  factor 
is  used  in  the  numerator  of  ilfj^(s)  while  the  same  Schuler  factor  is 
used  in  the  denominator  to  yield 


l_xm  1  xm 
bs  c 

(+) 


(4.45) 


Substituting  (4,45)  into  (4,44)  and  assuring  no  left— half  plane 
poles  at  the  roots  of  (+)  gives 


[II®  III'  »  lO(T>''s  4  i'')|  •  (ju  )  .  (,  4f, 

IT  ,.-|4j,.^ 

From  (4,46)  both  b  and  c  can  be  found. 

The  mean-squared  estimation  error  comparisons  are  accomplished  using 
the  basic  formula 


X  X 

NI 

Then  is  lound  by  siilistituti  ng  (4.45)  and  the  as  yet  undetermined 

i|f^(.s)  info  (d.SH).  Now  the  complete  i|r(s)  is  used  in  (4.42)  so  that 
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the  quadrature  in  (4.47)  may  be  mimcrlcall y  evaluated.  There  are  two 
distinct  parts  to  the  analytic  lortiiula  I  oi*  rms  position  estimation  error, 
one  from  the  Schuler  frequency  disturbances  and  one  from  the  low  fre¬ 
quency  disturbances  expressible  as 


) 


(4.48) 


These  are  found  from  substituting  (lr(s)  into  that  portion  of  (4.42) 
that  has  left-half  plane  poles,  after  (4.43)  is  satisfied.  Usinn  (4.42) 
modified  in  this  manner  as  the  into^rand  in  (4.47) 


—  i 

2«.1  J 


<1>^  (  s) 

S 

(+)(-) 
.t 


(4)(-) 


.j.(s)h'^,1,'^(h) 


ds 


(-) 


1 

2i',1 


_(s)uV''(.s) 

X  _ 

(-) 


(4.49) 


ds 


where  the  first  integral  represents 


,  AF  .2 

X 


X  (2nj) 


and  the  second 


) 


(-2H,))  . 
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Uslnti  ('1. 45)  ii\  the  ('a«ie.l»y  rcsi<liie  evaluation  ol  the  first  inteRral 


in  (4.49)  (jives 


,  AF  ^2 

<»«.  ' 

X  s 


2k(>) 


2o) 


(4.50) 


AK  2  2 

where  this  expression  is  equal  to  (Cr~  )  in  (4.18)  divided  by  w  .  This 

nV  ^  s 

is  a  standard  form  for  the  Schuler  mean-squared  estimation  error.  The 

form  for  the  low  frequency  counterpart  to  this  quantity,  however,  has 

no  such  standard  form  but  varies  with  the  number  and  location  of  low 

AF 

frequency  poles.  For  this  reason  the  quantity  (rr-  )  will  be  derived 

rSRx  t 

on  an  Individual  basis  for  the  measurement  combinations  specified  in 
Table  IV-1. 


Tl\e  external  velocity  measurement  may  be  included  as  a  measurement 
by  writ  I  HR  as 


(sXsfiR..,  - 

X  XX 

which  (jives  an  estimate  error  equation  of  the  form 


(4.51) 


SH  (s) 

X 


i|r(s)nr  -I  ,!f(.s)<:;  4 


4'(s)(j^ 

(4) 


l4<l>(s) 

S 


k»(s)1  yi 

(4)  -J 


(4.52) 


whe  re 

Ks)  ilr’(s)s  -  1  .  (4.53) 

WrltlnR  the  estimate  error  power  spc'ctral  density  for  the  position  error 
in  thf-  same  matrix  formal  ns  was  list'd  in  (4,2  )  yields 
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*8l», 


(a)  =  ^*(-8).  '^(-s)<[W2] 


♦  (a) 


*(s) 


2  2 
S  *li) 

a^(-) 


1 


s2  _  l<»(s)  \ 

(O  ' 


^  (4.54) 


where 


[W,]  = 


H4.^  (.s)h'^+R 

rr 


H«>-  (s) 

To 


(s) 

tv 


TT 


.  2  2,2 
(s  ) 
s 


(+)  (-)  .s^(4)(-)  ^ 


(4.55) 


Reasoning  identical  to  that  used  in  Sect.  A-2  yields  the  following  Schuler 
frequency  form  for  <l)(s): 


4.(s)  =  - 


(+) 

(+') 


(4.5f) 


Notice  that  the  low  frequency  behavior  of  (4.56)  gives 


'(s)|  ?•  -  1 

small  s 


(4.5'  ) 


This  results  in  highly  desirable  behavior  because  it  ininiitii/.es  the  con¬ 
tribution  of  the  1  +  'tifsl/s  term  in  (4.52).  This  term,  for  small  s 
and  1)(s)  ^  -  1,  might  be  very  large  even  for  reasonable  levels  of 
external  velocity  measurement  noi.se.  Therefore  it  will  be  assumed  that  the 
BVjj/s  term  in  (4,52)  is  large  at  low  frequency  making  l>(s)  -  -  f(4)/(+')] 
an  ideal  choice  at  both  low  and  Schuler  frequency. 

.Substituting  this  value  for  ''’(s)  into  (4.52)  yields 


W<,(s) 


ijr(s)lir  4  ^ll(s)p 


g 

X 


(I  ') 


k’fiv'" 

_ X 

(4  ') 


(4.58) 
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which  is  identical  to  (4.41)  except  that  It  is  replaced  with  k' . 
The  I’e  fore 


ilr  (s) 
s 


l)S  -f  c 
(+’ ) 


where  li  and  c  are  determined  from  tlie  condition 


rH4> 


IT 


— T  — T 
R](bs  +  c  ) 


«4g 

X 


(4.59) 


"^2 


4.1  U) 


The  estimate  error  PSD  is  expressible  as 


„r(-s)[H$^.^(s)n'"  ,  hV(s)  4 


(4.60) 


where 


>lf(-s)H$l  (s)  ^(s)h\’'(s) 

K  r 

X 


(+') 


(-’) 


2  2 
(- '  )  -  s  -  k '  s  4  , 


(4.61) 


AF 


NI 


and  n'~  and  a~  are  determined  as  previously  explained  Except  that 

^^x  X 

again  k'  is  used  in  place  of  k.  Tlie  resulting  Cauchy  residue  eval- 
.  AF  . 

nation  of  (CT--;,  ),  gives 

nit  S 
X 


.  AF  .2 

>s 

X 


B 


k’q 


2  '  2 

2k’-,,  2-. 

s  s 


(4.62) 


bl 


AF  2 

where  again  this  expression  is  related  to  (0'~„  )  in  (4.37)  by  the 
2 

lactor  CO  . 

s 
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In  the  next  two  sections  the  various  lorms  for  \|r  (s)  lor  the 

I 

cross-track  and  in-track  directions  will  be  found.  These  filter  forms 
will  be  constant  whether  or  not  the  external  velocity  measiii'oment  is 
used  as  a  measurement.  This  is  a  direct  result  of  choosint^  tfs)  as 
in  (4.56)  which  gives 


ilf*(s)  =  0  . 

small  B 


(4.63) 


Simply  stated  the  external  velocity  measurement  is  assumed  to  add  very 
little  low  frequency  information  and  therefore  does  not  change  the  low 
frequency  part  of  the  gradiometer  measurement  filters. 


B-2  Low  Frequency  Filters.  Cross— Track 

The  reason  for  the  low  frequency  filter  variation  from  cros.s-track 
to  in-track  is  solely  a  manifestation  of  the  assumed  low  frequency 

forms  for  ♦  (s)  and  (.s)  from  the  Heller  model.  Tl^e  low  frequency 

^y 

PSD  (power  spectral  density)  matrix  in  (3.11)  transformed  to  a  time 
variation  is 


2 

-k  s 
1 


C 

-k2S 


2k. 


-fik. 


C 


2k. 


-6k. 


-6k. 


24k.. 


-2k!|s 


-6k. 


24k. 


(4,64) 


A  B  C 

where  kj ,  k^ ,  and  k^  ar«  low  frequency  time  spectral  density  coeffi¬ 
cients  defined  as 


(4.65) 


and 


k  , 


•i  ' 


B 


J  ■  „(5-.i) 


C 
k  . 

.1 


,('>-.)) 


(4.66) 


(4.67) 


(4.68) 


Krotn  (4.64), 


•Ti  (s) 

K 


and 


s) 


are  of  degree  two  and  zero,  respectively. 


Thi.s  dilfercnce  ultimately  produce.s  low  frequency  filters  of  orders  two 
and  one.  Tliis  means  that  the  results  from  this  and  the  next  section  are 
applicable  to  any  other  gravity  perturbation  model  whose  low  frequency 
behavior  lor  4)  (s)  and  4*  (s)  matches  that  given  in  (4,64). 


The  cros.s-track  low  frequency  filter  derivation  begins  by  writing 
(4.43)  for  small  s.  Then 


r h/  .-(s)!!*^  h  K>^(s) 
rr 


"4 


(s) 


Sx. 


(4 . 69) 


small  s 


must  have  no  left-half  plane  poles.  To  satisfy  (4.69)  it  is  immediately 

T 

clear  that  tlie  denominator  of  ilf^(s)  must  be  the  left— half  plane  roots 
of 


lllip'  .j,(s)  H*  4  kI 

ff 


small  s 


(4.70) 


In  the  cross-track  direction  this  results  in  a  single  real  root  annotated 

s  -  .  (4.71) 
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To  assure  that  the  low  frequency  peak  in  the  estimate  error  sped  rum 

dies  out  at  the  upper  end  of  the  low  frequency  rantte,  the  nup'oratoi-  of 

ilf  (s)  must  attain  be  a  well  cljosen  constant  (s) .  This  is  equivalent  to 

satisfying  the  position  error  est'.mation  equivalent  of  the  convergence 

condition  in  (4.12),  at  low  frequency.  Therefore  \lf.(s)  lias  the  lor'ii 

X/ 


s  +  ^ 


(4.72) 


To  find  a  note  that  (4.69)  must  not  only  be  without  left-half  plane 
poles  but  to  assure  convergence 


-T 

.(s)h'^] 


”  ^  "large  s" 


(4.73) 


where  the  quotes  indicate  "large  s"  in  the  low  frequency  range.  Note 

also  that  (4.73)  is  applied  only  to  the  higher  powers  of  s.  It  may  be 

t  2  t 

assumed  that  ffi--^(s)  has  a  dominant  s  term  and  that  (s)  has 

a  dominant  s  term  when  estimating  position  error  cross-track.  To 

reduce  (4.73)  to  one  nontrivial  equation  using  these  facts,  write 


':'__t(s) 


2  t” 

s  .j.(s) 
TV 


s<t>-  (s) 

Fg 


(4.74) 


(4.73) 


where  the  primes  indicate  differentiation  with  respect  to  s.  Note  that 

(4.74)  and  (4,75)  only  hold  for  the  highest  powers  of  s.  Substituting 

(4.74)  and  (4,75)  into  (4,73)  yielils 
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(4.76) 


2 


IIO*' 

rf 


-T 

n 


1 


(«) 


L  ' 

^  Pk 

(ii  V 

s 


"larfre  s" 


Ki'oi.'i  (1.76)  one  I'oiiipotK'nt  oT  O  may  be  founc).  The  other  components  r'.ay 
be  iound  by  enforcinff  (4.69)  at  the  low  frequency  pole.  This  relation¬ 
ship  is 


^(s)  +  111 


— T 

rr 


=  0 


(4.77) 


Ucc-all  that  to  caUnilate 

evaluated.  Usiivu  the  Cavic 
jjives  no  standai'di/ed  form 


(0'~  )  this  last  inte};ral  in  (4.49)  must  be 
'  ^ 

by  residue  theorem  to  perform  this  quadrature 

AF 

as  was  the  case  for  (ff^  primarily  be- 


X 

cause  s  at  low  frequency.  Because  of  this 

and  the  nontrivial  nature  of  the  low  frequency  filter  constant  determin¬ 
ations,  detailcul  d('rivations  for  the  in— track  and  cross-track  low  fre¬ 
quency  filter  constants  arc  {{iven  in  App.  C.  These  derivations  apply 

to  measurement  combinations  listed  in  Table  IV-1  and  include  the  result- 

ak 

inu  detailed  formulas  for  . 


X 
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The  denominators  of  the  in— track  filters  are  aKain  found  from  (d.TO). 
For  the  assumed  second  decree  form  of  ^  (s),  however,  the  left-half 

plane  roots,  after  some  approximations  explained  in  App.  C,  turn  out 
to  be  a  complex  conjugate  pair  expressed  as  the  factor  (jC+).  The  exact 
location  of  these  poles  vary  with  the  exact  measurement  combination  used 
and  are  given  in  App.  C.  The  low  frequency  filters  can  be  written  as 


1 xm  1  xm 

t/s)  =  -  -  (-1.78) 

*  (ft) 

where  the  first-order  numerator  satisfies  convergence  at  low  frequency 
while  ensuring  sufficient  estimation  accuracy.  The  convergence  and  low 
frequency  pole  conditions  (Sect.  B-2)  are  used  to  find  ri  and  e.  These 
conditions  are  stated  below  for  the  in-track  direction 


(  ff  (/+) 


0  . 


(4.79) 


and 


.„(s)H^  +  m  (d^s  4  e^)|  ,  -  0  .  (4.80) 

— T  's:- roots 

TT  r  / .  ^ 

of  (f4) 

The  details  of  finding  d  and  e  are  again  peculiar  to  the  particular 
measurement  combination  <-onsldered. 


C  .  UOW  FREQUENCY  FIT.TEIl  KESUl.TS 
The  re.sults  from  the?  low  f  r<'<(U('in’,  y  I'iltci-  dt'ri  vat  ions  i  tt  App.  C 


are  given  in  Tables  IV-2  :ind  IV-:t. 


Tahlc  IV-2 

LOW  FRKQUKNCY  I'lLTKU  IlFSULTS 
(CUOSS-TRACK  MEASUREMENT  COMBINATIONS) 


These  mcasiircnient  combinations  ass\ime  the  state  vector  defined 
In  (A.l) 


Table  IV- 3 

LOW  FUF^UlvNCY  FILTER  RESULTS 

(in_track  measurement  comuinations) 


These  measurement  combi nal ions  assume  Liu?  state  vector  clel'lncd  In  (4.1) 
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Note  thill  to  nhtnin  the  totnl  analytic  position  error  variance 

formula  for  any  of  the  combinations  in  these  tables  (4,50)  or  (4.62) 

AP  o 

should  lie  added  to  the  appropriate  formula  for  (CT  y.  )  . 

OR 

X  i 


Cliaptcr  V 


REAL  TIME  KILTER  COVARIANCE  STUDY 


A.  COVARIANCE  COMPARISONS  OK  KAIAIAN  VS  APPROXIMATE 
KILTERS  KOR  THE  CROSS-TRACK  DIRECTION 


In  this  section  the  second— order  gravity  model  Is  used  in  a  standard 


Kalman  covariance  analysis.  A  comparison  is  then  made  with  the  covari¬ 
ance  results  obtained  from  the  previously  derived  approximate  filters. 
The  cross-track  direction  is  chosen  because  the  low  frequency  behavior 


The  cross— track  gravity  perturbation  power  spectral  density  from 
the  second-order  model  is 


<l>  ( s)  = 


(5.1) 


From  this  rational  expression  it  follows  that 


<1-;  (s) 

r  g 

xy  y 


X 


xy  xy 


Ms-'  -  nY 


)■ 


(5.2) 


(5.3) 


The  Kalman  state  space  representation  is  obtained  by  combining  (2.C), 
(2.8)  and  (3.3)  as  follows: 


X  K  X  G  w 


where 

X  state  vector 

F  state  transition  matrix 

G  process  noise  transition  matrix 

w  process  noise  vector  , 

The  ineasureincnt  is  asstimed  to  f)e 

z  =  H  X  4  e 


(5.5) 


where  T1  is  the  tadiomcter  error.  II  the  external  velocity  measurement 
is  Inclnded 


Tlie  first  measurement, 


here,  l.s  obviously 


which  can  be  written 
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(5,7) 


K  =  r  V 

y  xy  X 


which  comes  from  the  following  identity 


r  .  V  = 


r«  1 

rr  V  n 

X 

XX  X 

g 

r  V 

y 

xy  X 

g 

r  V 

L  zj 

Ixz  xj 

(5.8) 


V  =V  =0 

y  z 


Comparison  of  the  measurement  noise  in  (5.(>)  and  the  process  noise  in 
(5.4)  shows  that  the  two  are  correlated.  The  standard  Kalman  filter 
covariance  analysis  can  still  be  made  by  deriving  an  equivalent  system 
described  by  Bryson  and  Ho  [Ref,  16,  pp.  371  and  402-404]  where  the 
process  and  measurement  noises  are  not  correlated.  The  general  form 
for  the  equivalent  system  is 


X  =  F*X  +  G’w'  +  Dz 


(5,9) 


where 


D  = 

T  -1 

GT  R  (so  that  the  noise  in  (5,9)  is  imcorrelated  with 

g) 

(5.10) 

F'  = 

F  -  DH 

(5.11) 

G»  = 

r,  -D 

P  ^  ^ 

(5.12) 

w'  = 

[-:] 

(5,13) 

T8(t  - 

T)  =  E[i(t)  ;'^(T)] 

(5.14) 

The  specific  form  for  these  equations  using  the  unaltered  measurements 
given  in  (5.6)  is 
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Equation  (5.7)  shows  that  using  r  as  a  measurement  Is  equivalent 

X  y 

to  the  use  of  This  means  that  the  covariance  results  obtained  by 

slightly  modifying  the  CT— 1  formulas  in  Chapter  IV  can  be  compared 
with  the  covariance  results  from  the  foregoing  Kalman  configuration. 
Because  the  low  frequency  behavior  of  (s)  Is  different  for  the 

second-order  model  (sOM)  and  Heller  model^  (HM)  as  follows: 


.HM 


(s) 


=  2k„ 


small  s 


(5,20) 


SOM 

(=) 


small  s 


(5.21) 


the  low  frequency  root  from  (4.70)  is 


(5.22) 


The  value  of  a  found  by  substituting  (5.2)  and  (5.3)  into  (4.76)  is 


V 

a  =  .  (5.23) 

U) 

s 


The  analytic  low  frequency  variance  formula  is  altered  to  give 


q^a  -fi 

~  2 
3  V 

X  s 


(5.24) 


The  previous  results  for  the  determination  of  a,  b,  and  c  remain  unalt¬ 
ered  as  do  the  variance  formulas  except  that  the  spectra  In  (5.1),  (5.2), 
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and  (3,3)  are  used  iu  place  of  those  from  the  Heller  model. 


T’lere  arc  primarily  three  questions  to  be  answered  concerning  the 
real  time  approximate  filters:  (1)  How  optimal  are  the  approximate 
filters?  (2)  Over  what  range  of  carrier  vehicle  velocities  can  these 
filters  be  used?  (3)  What  measurement  combination  produces  the  most 
accurate  estimates  of  position  and  velocity  error?  Tables  V-1  through 
V-4  seek  to  answer  the  first  two  questions  by  comparing  the  Kalman  filter 
rms  estimate  errors  with  these  same  quantities  resulting  from  use  of 
the  approximate  filters.  The  assumed  second— order  model  and  inertial 
navigation  system  quantitative  characteristics,  taken  from  Heller  [10] 
are  as  follows; 


6V 


111 

y 


0.2  percent  of  V  -  band  limited  white  noise  with 
x 

a  two  second  correlation  time. 


(5.25) 


1,4466  X  10  ^  m^/sec^ 

7.5  X  10  ^  l/sec  , 


(5.26) 

(5.27) 


Because  the  Bell  and  Draper  gradlomoters  are  accurate  to  approximately 
1  Eotvos,  this  would  seem  a  realistic  value  to  use  and,  in  fact,  all  of 
the  results  in  this  research  assume  a  gradlometer  with  this  accuracy. 

In  addition,  the  gradlometer  is  assumed  to  have  a  10  sec  averaging  time. 
This  deteriiilnes  tlie  gradlometer  measurement  error  power  spectral  density 
necessary  for  computation  of  the  Kalman  filter. 

Looking  at  Tables  V-1  to  V-4  the  most  Important  results  appear  in 
the  rms  estimate  error  columns.  As  explained  earlier  the  dominance 
assumptions  made  in  the  filter  derivations  may  be  directly  checked  (which 
indirectly  chocks  optimality)  by  comparing  the  analytically  and  numerically 
obtained  rms  estimate  error  values.  The  last  two  columns  of  each  Table 
contain  these  numbers.  A  comparison  of  these  values  shows  that  nowhere 
do  they  differ  by  more  than  6%  and  routinely  differ  by  less  than  1%. 

A  direct  check  of  optimality  may  be  done  by  comparing  the  Kalman  filter 
estimate  error  results  in  column  2  with  the  numerically  Integrated  values 
in  the  last  column.  The  fact  that  these  values  also  agree  to  at  least  6% 
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and  routinely  to  less  than  1%  provides  direct  confirmation  of  approximate 
filter  optimality. 

Although  the  data  in  these  Tables  look  very  promising  it  must  be 
remembered  that  a  rational  second-order  model  was  used.  The  primary 
purpose  of  Tables  V-1  through  V-4  is  to  show  that  confirmation  of  the 
dominance  assumptions  made  in  the  approximate  filter  derivation  coincides 
with  a  near  optimal  filter.  From  this  point  on  then,  confirmation  of 
these  dominance  assumptions  will  be  considered  sufficient  proof  of  near 
optimality  since  direct  Kalman  comparisons  are  impossible  when  the  Heller 
Model  is  used  in  the  next  section. 

Table  V-I 

VELOCITY  ESTIMATION  COVARIANCE 
COMPARISON  (E/EVM)*:  EALMAM  VS  APPROXIMATE 

KF  B  Kalman  Filter,  AF  e  Analytic  Formula,  NI  =  Numerically  IntCKrated 


Carrier 

Velocity 

(a/aee ) 

Kalman 

Filter 

Reaults 

n 

’Sr 

7 

(mm/aec) 

Approximate  Filter 

AF 

•  Vgy 

y 

(m/scc)  <mm/soc) 

Reaults 

NI 

“^tv 

y 

(aa/aec) 

0.1079 

0.514 

0.1079 

0.1066 

JM 

1.079 

5.14 

1.079 

1.076 

10.79 

51.4 

10.78 

10.77 

MM 

103.38 

472.4 

103.43 

103.16 

*  Bxcludlnic  the  external  velocity  neaaurement  (EVM) 

Table  V-2 

VELOCITY  ERROR  ESTIMATION  COVARIANCE 
COMPARISON  (I/EVM)*;  EAIMAN  VS  APPROXIMATE 


Kalman 

Approxlnete 

Filter  Reaults  I 

Carrier 

Filter 

Velocity 

(m/aec) 

Results 

d^(mm/sec) 

y 

• 

(ai/sec) 

k' 

(l/sec) 

AF 

®6V 

(m/acc) 

NI 

*SV 

(e^/sec) 

O.Ooll 

0.514 

6.376"“^ 

0.0519 

0.0443 

fm 

0.5118 

5.14 

6. 376''“^ 

0.5196 

0.4807 

5.118 

51.4 

6.376''-‘ 

5.196 

4.88 

50,0ti 

472.4 

50.86 

47.71 

*I/EVM  «=  Inclwtlnfc  th«  vvloclty  iMasurcment 


—  5 


Table  V-3 

POSITION  ERROR  ESTIMATICW  COVARIANCE 
COMPARISON  (E/eVM);  KALMAN  VS  APPROXIMATE 
KF  =  Kaliiiaa  Kilter,  AF  =  Analytic  Formula,  NI  =  Numerically  Intesrated 


Carrier 

Velocity 

(ni/sec) 

Kalman 

Filter 

Results 

a  ..  (mm/ sec 

oV  ' 

y 

Approximate  Filter 

Results 

H 

NI 

(m) 

0,5144 

9,79 

-3.35^^ 

-2.5l'^^ 

3,35^^ 

9. 79 

9.79 

5.144 

17,4 

-2,5l'^^ 

3.35'^® 

17,44 

17.43 

51.44 

32 . 04 

-3.3^^ 

-2.51^^ 

3,  35'^^ 

32,17 

32.11 

514.4 

90.95 

03.074^® 

-2.12^'* 

3.35' 

99,  87 

97,82 

Table  V-4 

POSITION  ERROR  ESTIMATKW  COVARIANCE 
C(»»PARISON  (I/EVM);  KALMAN  VS  APPROXIMATE 


1 - 

Carrier 

Velocity 

(m/sec) 

Kalman 

Filter 

Result! 

KK 

'"or 

y 

(m) 

Approximate  Filter  Results 

b 

(m-sec) 

c 

(m) 

a 

(m-sec ) 

k' 

(l/sec  ) 

AF 

y 

(m) 

H 

0.5144 

9.23 

-2.046^® 

1 

« 

00 

2.646'^® 

6.376''^^ 

9.238 

9.231 

5.144 

16.41 

-2,040^® 

-1.687'''“^ 

2.646’^® 

6.376^“'* 

16,43 

16.42 

31,44 

29,31 

-2. 040^^ 

-1,687'^'* 

2.646'^^ 

\-4 

6.376' 

29.45 

29.39 

514.4 

59.98 

-2.470^® 

-1.389^'^ 

2.646^® 

v-4 

6,11' 

63,47 

61.61 
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SINGLE  GRADIOMETER  MEASUREMENT  ANALYTIC 
VS  NUMERICAL  COVARIANCE  COMPARISONS 


The  next  step  in  proving  the  accuracy  of  the  approximate  filters  is 
to  show  their  near  optimality  when  the  more  realistic  Heller  model  is  used. 
This  is  done  by  making  the  same  type  of  comparisons  as  in  Sect.  A  but  with 
the  following  changes:  (1)  The  Kalman  covariance  is  excluded;  (2)  Tie 
exact  formulas  for  Ch.  IV  are  used  (3)  The  IT-1  measurement  combination 
is  included.  The  Impetus  to  improve  the  Heller  model  came  from  studies 
done  in  this  section.  Therefore  all  the  comparisons  in  Tables  V-5  to 
V-12  were  done  assuming  the  BHM  (Baseline  Heller  Model).  In  addition 
these  data  also  assume  the  EVM  (external  velocity  measurement)  noise  nnd 
the  value  of  k  as  defined  in  Sect.  A. 

The  primary  focus  in  Tables  V— 5  through  V-12  is  the  agreement  of  the 
rms  estimate  error  values  in  the  last  two  columns  of  each  Table.  As  with 
the  preceding  Tables  in  this  Chapter  these  numbers  agree  to  within  6% 
for  both  velocity  and  position  error,  with  the  exception  of  the  lowest 
velocity  case,  in  all  of  the  velocity  error  Tables,  This  is  due  to  the 
fact  that  the  low  frequency  portion,  as  assumed,  of  the  estimate  error 
power  spectral  density  (PSD)  is  dominant.  The  exponential  decrease  in 
the  Heller  gravity  perturbation  spectrum  results  in  very  small  gravity 
perturbation  PSD  values  at  the  short  spatial  wavelengths  corresponding 
to  Schuler  frequency  at  low  speed  which  in  turn  causes  the  low  frequency 
dominance.  In  order  for  the  velocity  filters  to  be  meaningful  at  low 
velocities  some  form  of  low  frequency  compensation  would  have  to  be  in¬ 
cluded. 

By  comparing  the  data  which  includes  and  excludes  the  external  vel¬ 
ocity  measurement,  it  is  possible  to  determine  the  impact  of  using  this 
potential  data  source.  The  rms  estimate  error  values  in  fable  V-7  vs 
Table  V— 8  and  Table  V— 11  vs  Table  V— 12  show  that  the  external  velocity 
measurement  contributes  little  to  enhance  position  error  estimation. 
Conversely  Table  V-5  vs  Table  V-6  and  Table  V-9  vs  Table  V-10  demonstrate 
that  an  average  increase  in  accuracy  of  50%  is  posstl)lc  for  velocity 
error  estimation. 


Table  V-5 

IN-TRACK  VELOCITY  ERROR  ESTIMATION  COVARIANCE 
CaaPARISai  (E/EVM):  ANALYTIC  VS  NUMERICAL 
AF  =  Analytic  Formula,  NI  -  Numerically  Integrated 


Carrier  Velocity 

! 

AF 

NI 

a 

®0V 

X 

X 

(ni/sec  ) 

(m/sec > 

(mm/sec) 

(mm/sec) 

0.5144 

1.985^“^^ 

0.0126 

7.098 

5.144 

5.144 

1.078 

1.078 

51.44 

51.43 

10.78 

10.78 

514.4 

470.0 

103.16 

102.89 

Table  V-6 

Ilf-TRACK  velocity  ERROR  ESTIMATION  COVARIANCE 
COMPARISON  (I/EVM);  ANALYTIC  VS  NUMERICAL 


Carrier 

AF 

NI 

Velocity 

a 

k’ 

X 

X 

(ra/«ec  ) 

(m/ sec ) 

( 1 / sec ) 

(mm/sec) 

(mm/sec) 

0.5144 

3.96'^-'^ 

4.1'-^ 

7.098 

5.144 

5.144 

\-4 

6.376' 

0.5196 

0.493 

51.44 

51.43 

\-4 

6.  376' 

5.195 

4.93 

514.4 

470.0 

\-4 

6.095' 

50.8 

47.62 

Table  V-9 


CROSS-TRACK  VELOCITY  ERROR  ESTIMATION  COVARIANCE 
COMPARISON  (E/EVM);  ANALYTIC  VS  NUMERICAL 


Carrier 

Velocity 

(m/  sec  ) 

a 

(m/sec) 

AF 

"bV 

y 

(mm/ sec) 

NI 

^5v 

y 

(mm/sec) 

0,5144 

2.257V^^ 

0.0126 

4,088 

5.144 

5.143 

1,079 

1.077 

51,44 

5,142V 

10,78 

10.77 

514.4 

4 , 926^^ 

105,58 

105.53 

Table  V-10 

CROSS-TRACK  VELOCITY  ERROR  ESTIMATION  COVARIANCE 
COMPARIS(Mi  (I/EVM):  ANALYTIC  VS  NUMERICAL 


Carrier 

Velocity 

(m/sec) 

a 

(m/sec) 

k’ 

(l/sec) 

AF 

^6V 

y 

(mm/sec) 

NI 

^3v 

y 

(mm/sec ) 

0.5144 

2.575\-23 

4.5l\-^" 

\-4 

\-7 

1.382' 

4.089 

5.144 

5.143 

6.376' 

0.519 

0.493 

51,44 

5.142  V 

\-4 

6.376' 

5.195 

4.965 

514.4 

4 . 926\^ 

6.24^“"* 

51.39 

50.98 

The  IHM  (Improved  Heller  Model)  was  developed  after  noting  the 
excessively  long  characteristic  settling  times  associated  with  the  low 
frequency  filter  poles  for  some  of  the  cases  studied  in  this  section. 
To  best  see  this  problem,  Table  V-13  catalogs  these  settling  times  for 
the  CT-1  and  IT-1  (cros.s-track  and  in-track)  measurement  combinations. 

Table  13 


LOW  FHEQUENCY  POLE  CHAIlACTERISTIC  SETTLING  TIMES  FOR  THE  BHM* 


Note  the  settling  time  of  1.647  days  even  at  51,44  m/sec. 

The  most  desirable  way  to  deal  with  this  problem  is  to  update  the 
low  frequency  knowledge  of  the  gravity  field  thus  reducing  the  low  fre¬ 
quency  residual  in  the  Heller  gravity  perturbation  model.  This  is  pre¬ 
cisely  what  was  done  in  Ch.  II  to  obtain  the  IHM.  The  updated  character¬ 
istic  times  for  the  IHM  are  shown,  along  with  the  values  for  the  BHM  in 
Table  V-14, 
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Table  V-14 


LOW  FREQUENCY  POLE  CHARACTERISTIC  SETTLING  TIMES  FOR  THE  IHM 


Carrier 

Velocity 

(m/sec) 

CT-i 

Settling  Time  (days) 

IT-l 

Settling  Time  (days) 

(1/L) 

0.5144 

1157.74  (1647)* 

106.995  (176.9)* 

5.144 

36.612  (52.1) 

6.017  (9.95) 

51.44 

1.158  (1.647) 

0.338  (0.56) 

514.4 

0.037  (0.052) 

0.019  (0.0315) 

♦ 

numbers  in 

parentheses  indicate  BHM  values 

There  Is  some  Improvement  shown  In  Table  V-14  but  in  the  low  velocity 
cases  the  times  are  still  inordinately  long. 

C.  LOW  FREQUENCY  POLE  PLACEMENT  TECHNIQUE 

There  is  another  way  to  decrease  the  settling  times  but  it  is  less 
advantageous  because  acceptable  filter  performance  is  balanced  against 
estimation  accuracy.  This  simply  involves  moving  the  low  frequency  poles 
to  the  left  on  the  negative  real  axis  to  a  point  where  the  time  response 
is  acceptable.  These  new  pole  locations  are  nonoptlmal,  thus 
degrading  estimation  accuracy.  The  hope  is  that  estimation  accuracy 
sensitivity  to  low  frequency  pole  movement  is  small.  Table  V-15  shows 
this  relationship  for  a  velocity  of  51.44  m/sec  only.  The  results  for 
5.144  m/sec  and  below  indicate  that  the  extreme  root  relocation  necessary 
to  give  acceptable  low  frequency  behavior  produces  an  unacceptable 
degradation  in  estimation  accuracy.  At  very  high  carrier  velocities 
the  low  frequency  settling  times  are  already  on  the  order  of  hours  or 
even  minutes  and  do  not  require  severe  adjustment. 
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Note  that  for  the  cases  where  the  EVM  is  excluded  the  estimation  accuracy 
is  degraded  by  a  factor  of  2  to  3.  When  the  EVM  is  included  the  de¬ 
gradation  is  at  most  20%  and  may  be  as  small  as  a  few  percent.  Therefore 
pole  shifting  may  be  a  reasonable  way  to  hasten  filter  convergence 
particularly  if  the  EVM  is  included  as  a  filter  input. 

Another  way  to  see  the  effect  of  the  pole  shifts  is  to  look  at  plots 
of  the  estimate  error  PSDs  (power  spectral  densities)  for  approximately 
optimal  and  suboptlmal  pole  locations.  Figure  V-1  shows  these  PSDs 
for  the  IT-1  measurement  combination  excluding  the  EVM  assuming  a 
carrier  velocity  of  51.44  m/sec.  The  low  frequency  suboptlmal  PSD  is 
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FIG.  V-1  APPROX  I  MATE  I.Y  OPflMAL  VS  SUBOPTIMAL  POSITION  ERROR 

ESTIMATE  ERROR  POWER  SPECTRAL  DENSITY  COMPARISON. 


much  greater  than  the  same  portion  of  the  approximately  optimal  spectrum 
while  the  same  Schuler  peak  is  projliiced  in  both  curves. 


n.  MULTIPLE  GRADIOMKTER  MEASUREMENT  ANALYTIC 
COVARIANCE  COMPARISONS 


Up  to  this  point  only  the  CT— 1  and  IT— 1  measurement  combinations 
have  been  explored.  Do  the  other  combinations  provide  significant  im¬ 
provement  over  the  single,  pure  gradlometer  measurement  cases?  This 
question  is  answered  in  Tables  V-16  through  V-29,  The  covariance  data 
in  these  Tables  was  all  obtained  from  the  analytic  formulas  using  the 
IHM.  It  was  deemed  unnecessary  to  show  numerical  corroboration  since 
the  optimality  of  the  approximate  filter  technique  has  been  shown  in 
the  preceding  sections  of  this  chapter.  The  lowest  velocity  ease  is 
omitted  for  velocity  error  estimation  because  approximate  filters  are 
not  optimal  at  such  low  speeds,  as  explained  in  Sect.  B  of  this  chapter, 

Tlie  most  significant  information  occurs  in  the  last  columns  of  Tables 
V-16  through  V-29,  Two  types  of  comparisons  may  be  done  with  this  rms 
estimate  error  data, 

1,  The  different  measurement  combinations  may  be  examined  to  deter¬ 
mine  the  optimal  (lowest  rms  estimate  error)  combination.  This 
requires  that  the  estimate  error  values  in  the  1st,  3rd  and  5th 
rows  or  the  2nd,  4th,  and  6th  rows  be  compared  for  the  in-track 
Tables.  In  the  cross— track  Tables  the  numbers  in  the  1st  and 
3r<l  or  2n<l  and  4th  rows  should  be  examined.  This  process  clearly 
segregates  the  data  depending  upon  the  use  or  non-use  of  the 
external  velocity  measurement.  The  most  significant  conclusion 
is  that  the  IT-2  measurement  combination  is  inferior  for  In-track 
navigation  when  compared  to  the  IT-1  or  IT-3  combinations,  pai^ 
tlcularly  in  the  cases  where  the  external  velocity  measurement 
is  omitted.  Therefore  it  would  take  at  least  two  Bell  RGGs 
(IT-3  measurement  combination)  to  Improve  significantly  in— track 
navigation.  The  cros.s— track  data  shows  no  great  Improvement  when 
two  UGGs  (CT-2  measiirement  combination)  are  used  instead  of  only 
one  (CT-1  measurement  combination).  These  trends  confirm  the 
•ntuitlvely  obvious  fact  that  the  best  information  comes  from 
be  in-track  derivative  of  the  desired  gravity  perturbation  con^ 
ponent,  Pjjx  in-track,  and  cross-track.  Unfortunately 

one  UGG  will  supply  a  measurement  of  P^  but  it  takes  two  RGGs 
to  adequately  determine  Pj^j^.  ^ 
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The  improvement  in  estimation  accuracy  possible  throuRh  the  use 
of  the  external  velocity  measurement  can  be  found  by  comparing 
the  E/EVM  and  i/EVM  cases  for  each  measurement  combination. 

The  information  added  by  the  external  velocity  measurement  im¬ 
proves  velocity  error  estimation  by  an  average  of  50%  and  posi¬ 
tion  error  estimation  by  not  more  than  lO^’i. 


Table  V-16 

IX-TRACK(IT)  VELOCITY  ERROR  ANALYTIC  COVARIATSCE  C0UPARI30N 
FOR  =  5.144  a/aec 


Measurement 

Coablnation 

• 

(i»/sec) 

k* 

(l/sec) 

1 

AF 

'3v 

X 

(Msec) 

IT-1 

(E/EVM) 

SA44 

— 

1.079 

lT-1 

(I/EVN) 

5.144 

6.376'^"^ 

0.52 

IT-2 

(E/EVM) 

1.076^^ 

— 

91.95 

IT-2 

(I/EVM) 

1.076^' 

5.473'“"* 

4.814 

IT- 3 
(E/EVM) 

«.858l 

.-a.o  J 

— 

1.602 

IT- 3 
(I/EVM) 

1 

mil 

9.505^"'* 

0.634 

AF  X  analytic  formula 


Table  V-17 


IN-TKACK(IT)  VELOCITY  ERROIJ  ANALYTIC  COVARIANCE  COMPARISON 
FOR  V  “  51.44  m/scc 

X 


Measurement 

Combination 

a 

(m/sec) 

1 

k' 

(l/sec) 

af"^ 

''6V 

X 

(ram/sec ) 

IT-1 

(E/EVM) 

5,143^^ 

— 

10.  79 

IT-1 

(I/EVM) 

5. 143\* 

i 

v-4 

6.376' 

5.195 

IT-2 

(E/EVM) 

7.139 

— 

693.53 

IT- 2 
(l/EVM) 

7.139 

4.128'^”^ 

41.81 

IT-.'l 

(E/EVM) 

1 

— 

16.01 

IT- 3 
(I/EVM) 

1 

1 

[B 

V— 4 

9.503' 

e’.34 

AF  =  nnnlytic  formula 


Table  V-18 


IN-TRACK(IT)  VELOCITY  ERROR  ANALYTIC  COVARIANCE  COMPARISON 
FOR  V  -  514.4  m/sec 

X  ' 


Measurement 

Combination 

a 

(m/sec) 

k' 

(l/sec ) 

AF 

"’W 

X 

(mm/sec ) 

IT-1 

(E/EVM) 

3.423^^ 

— 

88.29 

IT-1 

(I/EVM) 

2 

3.423^ 

v-4 

5.202'“ 

46.93 

IT-2 

(E/EVM) 

5.216 

— 

149.4 

IT-2 

(l/EVM) 

« 

5.216 

\  —4 
8.842^ 

61.18 

IT- 3 
(E/EVM) 

1 

— 

109.11 

IT- 3 
(l/EVM) 

mm 

1 

6,45 

52.26 

Table  V-19 

CROSS-TRACK (CT)  VELOCITY  ERROR  ANALYTIC  COVARIANCE  COMPARISON 
FOB  V  =5.144  m/sec 

X  ' 


Measurement 

Combination 

a 

(m/sec) 

k' 

(l/sec) 

AF 

®Biv 

X 

(mm/ sec ) 

CT-1 

(E/EVM) 

5.143 

— 

1.079 

CT-1 

(I/EVM) 

5.143 

6.377'^’'^ 

0.520 

CT-2 

(E/EVM) 

— 

1.079 

CT-2 

(l/EVM) 

HI 

\“4 

6,379^ 

0.5197 

In  the  cro.sN-track  direction,  that  member  of  a  calculated  from 
Eq,  (4,14)  corresponding  to  the  measurement  of  Fy^  Is  an 
Imaginary  quantity.  This  indicates  that  the  numerator  of  the 
transfer  function  associated  with  Fy^  should  be  written  with 
a  first  degree  s  term  as 


Sp  s  1 

yz  's=jUg 


The  foregoing  footnote  applies  to  Tables  V-19  to  V-21. 
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Table  V-20 


CROSS-TRACK  VELOCITY  ERROR  ANALYTIC  COVARIANCE  COMPARISON 
FOR  V  -  51.44  m/sec 

X 


Measurement 

Combination 

a 

(m/sec) 

k' 

(l/sec) 

AF 

X 

(mm/sec) 

CT-1 

(e/evm) 

5.142^^ 

10.78 

CT-1 

(I/EVM) 

5.142'^^ 

6.375^“^ 

5.195 

CT-2 

(E/EVM) 

5.137^^ 

v-2* 

-J  1.656'' 

— 

10.78 

CT-2 

(l/EVM) 

mm 

1 

6,373''”'* 

5.194 

*  See  footnote  Table  V-19 


Table  V-21 


CROSS -TRACK  VELOCITY  ERROR  ANALYTIC  COVARIANCE  COMPARIS(»( 
F'OR  V  =  514,4  m/sec 

X  ' 


k» 

(l/sec) 


AF 

X 

(mm/sec) 


104.26 


6.161 


\-4 


51.07 


102.9 


.-4 


6.079 


50.73 


Measurement 

Combination 


CT-1 

(E/EVM) 


CT-1 

(I/EVM) 


CT-2 

(E/EVM) 


CT-2 

(I/EVM) 


a 

(m/sec) 


4,802 


,\2 


4.802 


\2 


4 . 674^^ 


-J2.079 


1.674^' 


-j2.079 


See  footnote  Table  V-19 
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Measuremen' 

Comblnatioi 


b 


c 


CT-1 

(E/EVM) 


(m-sec)  (m) 


-1.675^"^^ 

-6.291^“^'* 

-1.325''“^^ 

-2.993''"^'* 

CT-1 

(I/EVM) 


a 


(m-sec) 


k' 


(l/sec) 


AF 


(m) 


3.347''' 


10.58 


2.647V 


4.512 


X-15 


9,98 


3.347^° 


-1,314 


,\2 


9.48 


3.347 


.  2 


4.499 


,\-15 


9.48 


Table  V-2T 


CROSS-TRACK  POSITION  ERROR  ANALYTIC  COVARIANCE  COMPARISON 
FOR  =  5.144  m/sec 


Table  V-28 


CROSS-TRACK  POSITION  ERROR  ANALYTIC  COVARIANCE  COMPARISON 
FOR  V  =  51,44  m/sec 

X  ' 


Measuremenll 

Combination 

b 

(m-sec) 

— 

c 

(m) 

O' 

(ro-sec) 

k* 

(l/sec) 

AF  j 

y  1 

(m)  j 

1 

CT-1 

(E/EVM) 

-3.345^^ 

-2.508'^^ 

v7 

3,347' 

1 

34.58  j 

1 

CT-1 

(I/EVM) 

-2.646'^^ 

-1.686^^ 

2.646'^ 

S.-4 

6.375' 

31.78 

CT-2 

(E/EVM) 

-3.342^ 

5.172^“^^ 

-1.253''^ 

-1.336''“^ 

3.347'^ 

\5 

-1,313' 

31.2 

CT-2 

(I/EVM) 

f"  T 

-3.342'' 

-5.539''”^° 

r  \4'i 

-1.065' 

»1 

-1.336' 

_  - 

3.347'^ 

-1.313''^ 

\— 4 

6.373 

30.26 

J 

Table  V-29 


CROSS-TRACK  POSITION  ERROR  ANALYTIC  COVARIANCE  COMPARISON 
FOR  V  =  514.4  m/sec 

X  ' 


Measurement 

1 

AF 

Comblnatlot 

b 

c 

a 

k’ 

1 

y  1 

f 

1 

(tn-sec) 

(m)  ! 

1 

1 

(iB-sec) 

j 

(l/sec) 

1  (nt)  1 

!  1 

CT-1 

(E/EVM) 

1 

-3.124''® 

-2.187^^ 

3.347''® 

102.92  1 

1 

i 

CT-1 

(I/EVM) 

-2.5l''® 

-1.443'= 

\ -*4 

2.647' 

6.16l'‘~‘* 

j 

67.02 

CT-2 

(E/EVM) 

-3.04l'^* 

-5.16l''~^^ 

-1.140''‘'* 

-1.677^® 

3.347'® 

-4.152'®_ 

98.68 

CT-2 

(l/EVM) 

-3.041^® 

-5.84oV^^ 

-9.242'^^ 

-1.67?'^® 

3,347'‘* 

-4.152'® 

6.079''“^ 

67.26 

Chapter  VI 


GRAVITY  GRADIOMETER  SURVEYING  TO 
IMPROVE  THE  GRAVITY  PERTURBATION  MODEL 


The  gradiometer  may  be  used  as  a  survey  tool  divorced  from  any  direct 
INS  application.  Gradiometer  survey  data,  and  in  fact  any  survey  data 
which  provides  gravity  information,  may  be  used  to  improve  the  reference 
gravity  model  used  in  the  INS  (inertial  navigation  system)  and  consequently 
reduce  velocity  and  position  error. 

The  rms  velocity  and  position  error  can  b3  found  by  first  writing 
(4.4)  in  the  following  form 


m 


sg  skov”* 

IT)- "7^ 


(f.i) 


where  the  x  or  y  subscripts  have  been  omitted  as  the  equation  can  be 
used  in  either  the  IT  or  CT  directions.  The  velocity  and  position  error 
power  spectra  are  given  by 


/( s) 


k  q.. 


(  +  )(-)  (-H)(-) 


(P.2) 


and 


<v'=) 


s  O)  I,  s)  s  k  q.. 
K  S 


(^■)(-)  (+)(-) 


(6.3) 


Recalling  the  variance  formulas  in  (4.8)  and  (4.47)  the  position  and 
velocity  error  variance  can  oe  expressed  as 


;)R 


2 


2(<> 


+ 


_1_ 

IT.j 


ds 


(6.4) 


-HI- 


and 


s  4>^Asj 
(+)(-) 


ds  . 


(^.5) 


A(jain  tlie  intet;rnls  in  (6.4)  and  (6.5)  can  be  evaluated  either  numerically 
or  analytically.  Assuming  a  reduced  low  frequency  perturbation  content 
due  to  the  added  survey  Information,  the  analytic  integral  evaluations 
consider  only  the  Schuler  poles.  T)ie  validity  of  this  assumption  can 
again  be  tested  by  comparing  the  analytic  and  numerically  obtained  vari¬ 
ance  values.  Tl)e  analytic  formulas  are 


and 


s  s 


2k 


(6.6) 


(6.7) 


The  basic  theory  used  to  determine  the  updated  gravity  perturbation 
spectrum  lf<  taken  from  Breakwell  [l2].  It  is  presented  here  for 

completeness.  The  objective  is  to  estimate  some  scalar  function  -^l(x,y) 
related  to  the  surface  perturbation  potential  AU  by 

hj )  =  f((*>'’)uUv'‘J^;  O)  (6,8) 

where  (6.8)  is  written  in  the  Fourier  domain  and  lij  is  the  height  above 
sea  level  at  which  tlie  estimate  of  -T|  is  to  I)e  made.  The  measurement 
equation  may  Im’  expressed  as 


-82- 


(6,9) 


ttfXl  ntXl 

Z  (w®;  hg)  =  H  (w®)  AU(w®;  0)  +  i 

where  is  the  height  above  sea  level  at  which  the  survey  is  taken. 

Assume  that  the  estimate  of  at  height  h  is  given  by 


it  z(w®;  h^) 


(6.10) 


The  resulting  estimation  error  is 


h^)  =  h^)  h^)  (6.11) 

=  [  t'*'(  w®)h(  (I)®)]aU(  w®;0)+\|f^(  W®)e 


which  means  that  the  estimate  error  spectral  density  is  expressible  as 


X  Ch''''((Ii®)  t(w®)  -  +  t’'(-w®)R\|r(  io®)  . 


(6.12) 


In  this  situation  causality  is  not  a  factor  as  it  was  in  the  real  time 
filter  case.  Therefore  to  assure  the  minimum  mean  squared  estimation 
error,  given  by 


2 

CT~ 


h,  )d'‘)®dw®  , 
1  X  y 


(6.13) 


it  is  sufficient  to  minimize  (6,12)  with  respect  to  a(u®).  This  results 
in  the  statistically  best  weighting  of  the  measurements  for  estimating 
,  written  as 


ii  ill®)  =  [  n(  -“®)  >i’®y(  i"®  ;0)n  *  (  ■ 


‘n( !'2y( ■®:o)f("®)  .  (6.14) 
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Substituting  (6.14)  into  (6.12)  yields  the  following  compact  expression 
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1  +  3>A,r(^®;  o)h”^{w®)r"^h(-w®) 


(6.15) 


Notice  that  the  numerator  of  this  expression  represents  the  a  priori 
knowledge  of  an<l  the  second  term  in  the  denominator  represents  the 
information  added  by  the  survey.  To  obtain  the  required  one  dimensional 

g 

spectrum,  (6.15)  must  be  integrated  with  resjject  to  In  addition  the 

time  variation  spectrum,  required  in  (6.4)  through  (6.7),  is  found  simply 
by  dividing  the  one— dimensional  spacial  spectrum  by  the  IT  velocity. 

The  desired  spectrum  m.-iy  then  be  expressed  as 


Jo  i 


(6.16) 


Equation  (6.16)  may  be  u.se<l  to  determine 
(6.4)  and  (6.5)  11-*^  is  equal  to  g. 
the  updated  gravity  model  is  used  in  the 


the  updated  spectrum 

This,  in  fact,  is  the  case  when 

basic  navigation  loop. 


In  the  covariance  analyses  in  the  next  section,  two  survey  schemes 
are  used.  The  first  is  a  grouml  level  gradiometer  survey  ((SS)  where  the 
measurement  is  assumed  to  lie  the  vertical  perturbation  gradient,  r 
Recall  that  at  least  two  Bell  RGGs  are  required  to  obtain  this  quantity. 

In  all  the  a  priori  stmlles  the  BHM  is  taken  as  the  gravity  perturbation 
model  to  be  updated.  Therefore  the  surface  potential  fluctuations, 
from  (5.4),  are  descrllx?d  as 
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The  functional  relationship  between  (6,17)  and  r 


is  found  to  l>e 
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(6.18) 


and 


H(u^)H(-a)^) 
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The  second  type  of  survey  used  is  the  proposed  low-low  GRAVSAT  con¬ 
figuration  where  two  low  earth  satellites  follow  each  other  in  circular 
polar  orbits  and  are  connected  by  a  very  accurate  Doppler  rate  indicator 
The  relative  velocity  changes  between  the  two  satellites  indicate  gravi¬ 
tational  field  perturbations.  The  functional  relationship  between  the 
Doppler  data  and  the  potential  perturbations  is 
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(6.20) 


where  V  is  the  unperturbed  circular  velocity  and  A  is  the  satellite 
ir 

distance  spacing.  Assuming  that  the  perturbations  of  interest  have 
wavelengths  substantially  greater  than  /'  then  H((o''’)  from  (6.20) 
may  be  approximated  as 
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where  •©■  denotes  the  direction  of  u)'  .  Tlien 
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(6.22) 
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It  is  possible  to  modily  (6.22)  to  compensate  for  the  rotation  of  the  earth 
under  the  satellite  ('round  track  as  explained  by  Breakwell  [Ref, 

12,  p.  Tliis  correction  is  necessary  since  the  north-south  and 

south-north  satellite  ({round  tracks  are  not  parallel.  These  tracks  diverge 
at  an  angle  having  a  value  of  approximately  l/l6  radian  at  low 
latitude  .  To  compensate  for  this  effect  (6.22)  should  be  replaced  by 
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(•>  A  [cos"  '(e!--e)+cos  ]  (6.23) 
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or  simplifying  (6,23)  using  trigonometric  identities 
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(6.24) 
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For  both  the  gradlometer  and  satellite  surveys  the  measurement  error  is 
given  by  the  following  formula  in  the  scalar  case 


N/AS 


(6.25) 


where  is  the  mean-squared  measurement  error,  N  is  the  number  of 

survey  measurements,  and  AS  Is  the  area  surveyed. 

The  analysis  in  the  following  section  investigates  both  the  IT  and 
CT  directions.  The  transfer  functions  relating  g  and  g  to  the  BHM 
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perturbation  potential  are  respectively 
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With  these  relationships  the  error  spectra  from  (6,16)  and  subse¬ 
quently  the  INS  position  and  velocity  error  may  be  determined  for  the 
following  four  separate  direction/survey  options, 

IT/GS  (in-track/sradiometer  survey).  This  combination  Indicates  al¬ 
teration  of  the  spectrum  using  gradiometer  survey  data.  For  this 

case  (6.16)  becomes 
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(6,28) 


1T/SS  (in-track/satellite  survey).  This  again  indicates  alteration 
of  the  residual  but  this  time  with  a  low-low  GUAVSAT  satellite  survey. 

The  updated  g^  spectrum  may  be  written  as 
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where 


C(w®) 


N/AS 


a 


2 

Doppler 


•I- 


i 


-2 


e 


I 


(6.30) 


CT/GS  (cross-track/gradiometcr  survey).  This  again  Indicates  altera¬ 
tion  of  the 

survey.  The  updated  g 


gy  residual  but  this  time  with  a  low-low  GRAVSAT  satellite 
spectrum  may  be  written  as 
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CT/SS  (cross— track/satellite  survey).  Here  the  spectrum  form¬ 

ula  updated  by  satellite  data  appears  as 
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The  first  step  In  doing  the  covariance  comparisons  In  the  next 
chapter  Is  to  numerically  carry  out  the  quadratures  in  (6.28)  through 
(6.32)  to  obtain  the  updated  gravity  perturbation  spectra.  These  spectra 
are  subsequently  used  In  the  numerical  evaluation  of  (6.4)  and  (6.5) 
to  bo  compared  with  the  position  and  velocity  errors  obtained  from  (6,6) 
and  (6.7)  wliere  the  Schuler  values  from  the  updated  spectra  are  used  In 
the  analytic  formulas. 
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Chapter  VII 


COVARIANCF  COMPARISON  FOR  REAL  TIME  VS 
SURVEY  USE  OF  A  GRADIOMETER 


The  objective  In  this  Chapter  will  be  to  compare  the  real  time 
estimation  errors  in  Sect.  B  of  Ch,  V  with  the  velocity  and  position 
errors  resulting  from  the  use  of  an  updated  gravity  model.  To  obtain  a 
valid  comparison  the  survey  parameters  must  be  chosen  Judiciously, 
particularly  in  the  case  of  the  gradiometer  survey.  The  survey  vehicle 
velocity  and  gradiometer  averaging  time  are  the  principal  parameters. 
Assuming  that  the  latter  is  fixed  it  is  important  that  the  gradiometer 
survey  and  carrier  vehicles  have  the  same  speeds  to  assure  equal  wavelength 
resolution  in  the  real  time  and  survey  applications.  Actually,  the  im¬ 
portant  parameter  in  the  survey  situation  is  the  number  of  survey  measure¬ 
ments  per  unit  area  (N/AS),  If  the  survey  track  spacing  is  set  at  a 
distance  equal  to  the  minimum  wavelength  information  desired  (10  km), 
and  with  a  fixed  gradiometer  averaging  time,  the  survey  velocity  completely 
determines  the  measurements  per  unit  area  and  the  time  to  survey  a  speci¬ 
fied  area.  Then  the  claim  of  equal  survey  and  carrier  speeds  ensuring 
equal  resolution  is  contingent  on  the  track  spacing  chosen.  From  a 

practical  standpoint,  the  survey  times  for  a  given  area  are  of  interest. 

12 

Table  VII-1  gives  the  survey  times  for  an  area  of  1  X  10  square 
meters,  about  l/l5  of  the  continental  United  States. 


Table  VII-1 

SURVEY  TIMES  FOR  1  X  10^^  SQUARE  METERS 


Survey  Velocity 

Survey  Time 

(m/sec) 

(days) 

0.5144 

2250 

5.144 

225 

51.44 

22.5 

514.4 

2.25 
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Note  that  the  lower  velocities  take  excessive  survey  times,  Gradiometer 
survey  height  also  Influences  the  resolution  question.  Thus  to  maintain 
equality  the  ^radiometer  survey  and  carrier  vehicles  will  also  have  equal 
assumed  altitudes.  With  these  facts  in  mind  the  gradiometer  survey  para¬ 
meters  arc  ns  follows; 

survey  height:  zero  km 
survey  velocity:  carrier  velocity 
gradiometer  averaging  time;  10  sec 
track  spacing:  10  km 
gradiometer  accuracy:  1  E5tvos  , 

The  GUAVSAT  mission  parameters  are  more  interdependent  and  conse¬ 
quently  less  flexible  than  the  gradiometer  survey  quantities.  Nevertheless, 
the  GUAVSAT  survey  provides  an  interesting  alternative  to  the  gradiometer 
survey  with  very  real  possibilities.  Nominal  parameters  for  a  low-low 
GUAVSAT  mission  might  be 

satellite  altitude:  180  km 

circular  velocity:  7.796  km/sec 

-6  , 

Doppler  accuracy:  1  x  10  m/sec 
satellite  separation:  3  deg  (341.87  km) 
length  of  mission:  6  months 
time  between  measurements:  8  sec 

Assuming  these  survey  parameters,  Tables  VIl-2  and  VH-3  compare  the 
real  time  data  from  Tables  V— 5  through  V-12  with  the  variances  obtained 
from  (6.6)  and  (6.7).  In  order  to  standardize  the  data  the  optimal  value 
of  k  is  used  for  both  the  real  time  and  survey  applications.  In  the 
real  time  case  the  optimal  k  is  given  in  (4.31)  (x  or  y  subscript  as 
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Table  VI 1-2 

REAL  TIME  VS  SURVEY  VELOCITY  ERROR  COVARIANCE  COMPARISON 

EVM  =  External  Velocity  Measurement 
GS  =  Gradiometer  Survey 
SS  =  Satellite  Survey 


Measurement 

RMS  Velocity 
or  Velocity 
Estimate 

Combination 

Error 

(mm/sec ) 

AF 

4.1-^ 

0,5196 

5.195 

50,8 

IT-1 

(I/EVM) 

NI 

X 

7.098 

0.493 

4,93 

47,62 

AF 

4.095''"^ 

0,456 

2.574 

14.46 

it/gs 

X 

NI 

a 

fiV 

X 

0.0177 

0.422 

2,541 

14.32 

AF 

4.095'^"^ 

8.295 

16.94 

10.72 

IT/SS 

NI 

6.923 

5.95 

23.14 

11.54 

AF 

1.382''“^ 

0.519 

5.195 

51.39 

CT-1 

y 

(I/EVM) 

NI 

V, 

4.089 

0.4934 

4.965 

50.98 

AF 

‘"nv 

1.382'^“^ 

0.4022 

2.544 

14.46 

ct/gs 

y 

NI 

e.sse'^"^ 

0.386 

14.32 

AF 

1.382^'“^ 

5.18 

30.  37 

57. 18 

CT/SS 

y 

NI 

3.99 

4.46 

1 

22.77 

56.54 

Note:  Data  assumes  IMIM  (naseline  Ili'ller  Model)  prior  to  survey  update 
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Table  VI I- 3 


REAL  TIME  ys  SURVEY  POSITION  ERRQR  COVARIAHCf:  COMPARISON 

EVM  r  Ex ti'rnal  V»‘loc^Ly  MeasiifoiiicnI. 

GS  =  G 1- ad  i  OHIO  ter  Survey 
SS  =  Satellite  Survey 


RMS  Pqs  i  1 1  on 

CARR] 

E  R  V  E 

(m/sec) 

L  0  C  I  T  ' 

y 

or  Position 

Measurement 

Estlmaie 

Combination 

Error 

<m> 

Q.SH-i 

5.144 

51.44 

514.4 

AF 

3,79 

8.93§ 

21.47 

56,71 

IT-1 

(l/EVM) 

3,73 

3.71 

20.58 

58,76 

AF 

0.367 

2.076 

11.67 

IT/GS 

^^R 

(I.2I9 

0.556 

2.297 

11.72 

X 

AF 

-  .  V— 7 

— 

3.03' 

6.6^1 

13.06 

8.646 

IT/SS 

X 

NI 

"^fiR 

108.8 

11.55 

8.46 

X 

AF 

n.9 

CT  ^ 

m 

21,17 

37.83 

76.42 

c'r- 1 

y 

(r/EVM) 

NI 

11,89 

21.12 

37.54 

74.77 

y 

AF 

^R 

1,1 15^”^ 

0. 3244 

2.05 

11.67 

ct/gs 

y 

NI 

*^?iR 

0.113 

0.468 

2.23 

11.75 

y 

AF 

".R 

y 

' 

i.iisV"^ 

4.18 

24.47 

46.12 

ct/ss 

NI 

'"fiR 

110.82 

34.29 

25.94 

46.23 

y 

Note:  Tills  data  assumes  BHM  prior  to  survey  update. 
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applicable)  which  is  attained  in  practice  by  simply  lncludi>K  the 
external  velocity  measurement  in  the  IT-1  or  CT-1  combinations.  For 


survey  case, 


used  in  the  basic  navigation  loop  is  cliosen  iisins 


except  that  the  update  of  is  done  a  priori  so  that 


the 

(4.31) 
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where  the  value  on  the  right  hand  side  of  (7.1)  comes  from  the  updated 
gravity  perturbation  spectrum. 

From  Tables  VII-2  and  VI 1-3  three  things  are  immediately  apparent. 
First,  for  the  low  velocity  survey  cases,  agreement  between  the  analytic 
and  numerical  standard  deviations  is  poor  for  velocity  and  position 
error.  This  is  not  surprising  when  it  is  remembered  that  neither  the 
analytic  velocity  nor  position  results  take  into  account  any  low  fre¬ 
quency  contribution.  The  existence  of  a  dominant  low  frequency  contribu¬ 
tion  produces  the  observed  discrepancies.  The  second  obvious  fact 
that  the  gradlometer  survey  always  produces  more  accurate  results  than 
the  real-time  gradlometer.  Reflecting  upon  the  situation,  this  too  is 
is  an  expected  result  since  the  gradlometer  survey  gathers  information 
along  many  tracks  while  the  real-time  gradlometer  has  knowledge  only 
along  the  present  vehicle  track.  The  third  and  most  difficult  trend 
to  explain  is  the  GRAVSAT  vs  gradlometer  survey  behavior.  In  general, 

gradlometer  survey  gives  better  resvilts  than  GRAVSAT.  The  high  fre¬ 
quency  gravitational  potential  is  attenuated  with  increasing  altitude. 

This  makes  the  ground- level ,  1  Eotvos  gradlometer  more  sensitive  to  short 
spatial  wavelengths  (characteristic  of  Schuler  frequency  at  low  velocity) 
than  the  lO”®  m/sec  Doppler  at  satellite  altitude.  However,  a  second 
consideration  working  against  this  trend  is  that  as  the  carrier  velocity 
increases,  so  too  does  the  spatial  wavelength  corresponding  to  Schuler 
frequency.  It  .just  so  happens  that  for  the  highest  velocity  IT  case  the 
Schuler  spatial  wavelength  is  so  long  that  GRAVSAT  lias  a  superior  signal- 
to-noise  ratio  when  compared  with  that  of  the  gradlometer  at  ground  level. 
Had  speeds  above  514.4  m/sec  lieen  considered,  trends  in  the  data  suggest 
this  same  type  behavior  might  appear  in  the  CT  direction  also. 
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Chapter  VIII 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

_A.  CONCLUSIONS 

The  conclusions  resultinfi  from  this  work  can  be  classified  in  two 
areas:  those  resulting  from  the  real  time  work,  and  those  from  the 

survey  application  of  the  ^radiometer  to  inertial  naviKation, 

A-1  Real  Time  Gradiometer  Application 

It  has  been  shown  that  the  modified  Wiener  filter  analysis  produces 
low-order  filters  for  inertial  navipcation  system  velocity  and  position 
error  estimation. 

Analysis  of  the  accuracy  of  these  filters  Is  accomplished  in  several 
ways.  If  Kasper’s  rational  gravity  model  is  used,  covariance  compari¬ 
sons  with  the  Kalman  filter  show  agreement  in  the  rms  position  and 
velocity  error  estimate  error  typically  to  within  1%.  An  indirect 
method  of  proof  is  needed  for  the  more  accurate  Heller  gravity  model. 

This  is  done  by  a  comparison  of  analytically  and  numerically  obtained 
estimate  error  results.  Anreement  of  the  analytically  and  numerically 
obtained  estimate  error  values  is  always  within  lO'"',  on  the  averat^e 
about  5%,  and  in  some  cases  below  1%.  T)iis  assumes  that  for  velocity 
estimation  the  carrier  vehicle  has  a  velocity  at  or  above  5,144  m/sec 
(10  knots).  Therefore  for  the  single  ^radiometer  measurement  cases  used 
in  the  above  covariance  comparisons  the  approximate  low  order  filters 
are  very  nearly  optimal. 

The  approximate  nature  of  the  derived  filters  lends  to  the  carrier 
vehicle  velocity  limitation  mentioned  in  the  previous  para^rraph.  Fully 
stated  this  limitation  requires  that  the  position  and  velocity  filters 
be  used  only  on  carrier  vehicles  traveling  at  or  above  5.144  m/sec  (for 
the  assumed  levels  of  external  velocity  measurement  noise  and  ^radiometer 
accuracy)  to  avoid  at  least  two  problems  that  arise.  These  problem.s  are: 
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(1)  tlie  volocity  I'ilters  contain  Insul'f icient  lo\  frequency  gravity 
perturbation  compensation  witicli  makes  them  inaccurate  at  low  velocity 
vvliere  there  is  a  large  low  frequency  gravity  perturbation  contribution; 

(2)  The  position  error  filters  do  compensate  low  frequency  gravity 
disturbances  whicli  unfortunately  produces  excessively  long  low  frequency 
pole  characteristic  settling  times  in  the  low  velocity  cases.  In  spite 
of  a  pole  placement  technique  which  is  marginally  successful  for  mid¬ 
range  velocities,  and  the  improvement  of  the  Heller  model  to  reduce 

the  low  fi-equency  portion  of  the  gravity  potential  perturbation  spectrum, 
the  under  5.144  m/sec  velocity  range  is  unsuitable  for  approximate 
filter  application  at  the  present  time.  Lower  velocities  may  be  more 
suitable  as  future  gradlometer  and  GHAVSAT  surveys  provide  greater 
knowledge  of  the  earth's  gravity  field. 

The  ability  of  the  Wiener  technique  to  produce  approximate  filters 
for  multiple  measurement  situations  makes  it  possible  to  compare  the 
advantages  of  using  two  gradiometers  as  compared  to  only  one.  The 
effect  of  using  the  external  velocity  measurement  as  a  filter  input 
may  also  Ije  explored.  In  the  cros.s-track  direction  using  one  gradlometer 
(CT-1  measurement  combination)  as  opposed  to  two  (CT-2  measurement  com¬ 
bination)  produces  a  degradation  in  estimation  accuracy  of  not  more  than 
2''n  for  velocity  error  estimation  and  11%  for  position  error  estimation. 
However,  using  one  gradlometer  in-track  (IT-2  measurement  combination) 
as  opposed  to  two  gradiometers  (IT-3  measurement  combination)  yields  a 
degradation  in  velocity  error  estimation  accuracy  of  from  0.5  to  1.5 
orders  of  magnitude  and  in  position  error  estimation  accuracy  of  from 
1  to  2  orders  of  magnitude.  Clearly,  if  significant  Improvement  over 
the  a  priori  position  and  velocity  error  values  is  desired,  two  gradi¬ 
ometers  are  required  In-track  while  one  will  suffice  cross-track.  Addi¬ 
tional  estimation  accuracy  may  also  be  obtained  for  a  given  gradlometer 
measurement  combination  by  simply  Including  the  external  velocity 
mea.surement  as  a  filter  input.  In  most  cases  of  velocity  error  estima¬ 
tion  this  results  in  an  improvement  of  about  50%.  The  improvements  in 
position  error  estimation  accuracy  are  much  more  modest,  being  generally 
on  the  order  of  a  few  percent,  due  to  the  increased  effect  of  low 


Of  the  gradiometer  and  GKAVSAT  surveys  considered  In  this  research 
the  gradlometer  option  Is  more  accurate  by  at  least  a  factor  of  four, 
except  for  the  highest  velocity  in-track  case.  As  expected,  the  data 
show  that  the  gradlometer  survey  always  produces  more  accurate  results 
than  does  the  real  time  gradlometer  scheme.  However,  it  may  not  be 
possible  to  obtain  gradlometer  surveys  over  all  potential  carrier 
vehicle  routes.  Even  the  GRAVSAT  survey  with  worldwide  coverage  showed 
resolution  problems  at  spatial  wavelengths  greater  than  its  altitude 
caused  by  high  frequency  attenuation  of  the  gravity  potential.  There¬ 
fore,  in  spite  of  the  excellent  INS  accuracy  attainable  with  the  grad¬ 
lometer  survey,  the  most  likely  scenario  for  improving  INS  accuracy  in^ 
volves  immediate  real-time  use  of  the  gradlometer  with  a  gradual  improve¬ 
ment  of  the  gravity  model  from  data  collected  in  surveys.  These  data 
could  come  from  any  source  supplying  gravity  information,  even  the 
stored  real-time  gradlometer  outputs.  One  final  question  is  the  ability 
to  store  aboard  a  carrier  vehicle  the  large  number  of  coefficients 
necessary  to  describe  a  world— wide  gravity  model  with  a  desired  resolu¬ 
tion  of  approximately  10  km.  One  solution  to  this  storage  problem 
might  be  to  use  an  alternative  gravity  field  description  consisting  of 
strategically  placed  mass  concentrations  in  the  pro,iected  area  of 
vehicle  operation.  This  makes  it  unnecessary  to  store  a  world-wide  model 
thus  alleviating  the  storage  problem  while  making  the  actual  computa¬ 
tion  of  gravity  simpler. 

RECOMMENDATIONS  I  OR  FUTURE  RESEARCH 

The  approximate  filters  used  to  estimate  position  and  velocity 
error  from  the  real  time  gradiometer/lNS  combination  were  derived  by 
assuming  "white"  gradlometer  an<l  external  velocily  measiiremeiil  noises. 
Real  gradlometer  noise  is  colored.  Data  ali'eady  obtained  from  the 
gradlometer  manufacturers  and  proce.ssed  by  the  Analytic  Science  Corp. 
could  be  used,  along  with  Navy  te.st  data,  when  It  Incomes  available, 
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to  develop  more  discrltiiinatlht'  I'ilterst  l-’lttlnK  the  assumfed  external 
velocity  measurement  ndlse  to  a  particular  carrier  vehlfcle*s  measure¬ 
ment  equipment  would  also  lielp  to  develbp  i)etter  filters. 

The  metliorl  used  in  this  research  to  analyze  the  error  spectrum 
of  the  uravity  field  Improved  by  sUrvfey  data  provides  several  options 
for  further  study.  This  etrOr  spectrilm  analysis  method  could  be  used 
to  optimize  ^radiometer  survey  parameters.  Additionally  it  could  be 
used  to  compare  different  measurement  schemes  as  was  done  with  the  low- 
low  GRAVSAT  and  ground  leVel  ^radiometer  surveys  in  this  research.  Two 
other  survey  questions  that  mlf;ht  be  exhlored  are:  (1)  The  possible 
problems  and  merits  associated  with  storing  real  time  gradiometer/lNS 
data  for  post  processing  and  gravity  model  improvement;  (2)  The  advan¬ 
tages  of  flying  a  gradiometer  on  board  the  QtlAVsAf  spacecraft  as  an 
additional  data  source; 
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APPENDIX  A 


POISSON’S  formula  for  upward  continuation 

OF  THE  EARTH'S  PaPENTIAL 


To  prove  Poisson's  formula  It  is  necessary  to  relate  the  followln}; 


spherical  harmonic  forms  for  the  earth's  potential; 


(A.l) 


®  i  R,^  ^1 

t  *’ 

f|)  m=-£' 


(A. 2) 


where 


-  Vim**'”  «>• 


for  m  S  0 


(A.  3) 


(A. 4) 


Note  from  (A.l)  that 


»  i 


^  “<V  t  ^  iS  m5 


1  1 


(A. 5) 


Where  the  Integral  specified  is  over  the  entire  earth's  surface.  Assuming 

that  the  Y.  's  are  not  only  orthogonal  over  the  sphere  but  suitably  (fully) 
Jtm 

normalized  [17]  then 


(A.P) 
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SvibstitatinK  (A.e)  into  (AiS)  yields 

^  u(R,^  e,  »)?*^  (e,  ♦)«  =  ‘ 

There  Cofc 

^  ^  "‘V  »<  «>  *  • 

5?ubstltutinK  (A. 8)  Into  (A, 2) 

U(r,  B,  ♦)  -  I  ^  *’ 

X  <t>*)  u(Rp»  e'*4)*)dn‘  * 

Making:  us<‘  o(  ihv  idc^ntity 

£. 

Pg(cos  0  =  ^ 

where  ,]r  Is  the  onfjle  between  (I,  n,  <1))  and  (1,  O’,  (A. 9) 

be  written  as 


u(r,e,^)  =  If  p  (!i>)^'N2M)P/cds  Ou(R,^e‘>4>‘)«J0‘  . 


Since  i!f  is  the  an^le  betweert  R^  and  r  theh 


1«,.  ‘  ^1  "  I'-'*  t  ^  "T  ‘ 


UsirtK  the  neneratinK  liihetiort  foe  •'^<1*, )  expressed  aS 


[I  H.  =  3C  ^'"p  (I) 


and  substltutlrif?  ih  the  quantities 


(A.  7) 


(A  .  8) 


(A. 9) 


(A*10) 

cah 


(A. 11) 


(A. 12) 


(A. 13) 
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5  =  ana  §  =  cos  v 


then  from  (A. 12) 


[1  +  I*  -  2||«  ]■ 


Substituting  (A.  14)  and  (A.  15)  Into  (**.13) 


5  =  -  2R^r  cos  ^  E  (-jr’  )  Pg(cos  4^)  . 

Taking  the  negative  of  the  partial  with  respect  to  r  of  (A, 16) 


2rRj  «•  -  cos  'I')  ^  ' 


Subtracting  (A. 16)  from  (A. 17)  yields 


2rR  (r-R_cos  i|f)-R,D 

(B  (1)  ^1) 


R  R  .'-l-ll 


P^Uos\|f) 


or 


/  2  2. 
"(b''-  -  V 


% 

(2^l)(— )  Pg(cos  \ti)  . 


Putting  (A. 19)  into  (A. 11)  gives 


U(r,  0,  <t0  = 


/  2  2^ 


4jc 


U(R|/5S1>‘) 

(r^i-R^-2T.  R,/^2 


(A. 14) 

(A. 15) 

(A. 16) 

yields 

(A, 17) 

(A. 18) 

(A. 19) 

(A. 20) 


This  equation  applied  to  the  perturbation  potential  Is  the  desired  result. 


Appendix  B 


SPHERICAL  HARMONIC  AND  HELLER  MODEL  RMS  POTENTIAL 
COEFFICIENT  VARIATION  EXPRESSIONS 


The  GEM  10  [lA'l  data  is  given  in  terms  of  C  's  and  S^  's  from 

•'  jim  ?.m 


the  following  model 


“  ^  R| 

U(r,d,<t))  =  ^  S  m<J)+C^_^sin  nKt>]  .  (B .  1 ) 


£=0 


Therefore  the  equation  lor  (Eq.  3.15)  is  best  derived  using  the 

following  spherical  harmonic  form  first  defined  in  (A, 2) 


U(r,  0,  4>)  =  ^ 


£  R,  ^-1 


(B.2) 


is  is  first  necessary  to  establish  the  relationship  between  C  S 

fm  fm 

and  J  .  To  do  this  consider  the  plus  and  minus  m  pairs  in  (B.2) 
fm 

for  a  given  X 


R^  i+l( 

^  K._P 


«0  " 


{j  ,  K,  P  ,  (sin  O)  [  cos  m<J>  +  i  sin  mij)  J  i-  j  ,,  >.K  „ 

I  tm  m  ivm  i,(-m)  .(i( -n 


(B.  3) 


X  (sin  0)[cos  mi))  -  i  sin  m'l>  |  >  , 


Assuming  that 


P.;^(sinf*)  =  Pg^_^)(sin  0) 


(B.4) 


K  , 

i.m 


(-1) 


(2  •-  f,  )■ 
o  m 


(B.r.) 


-io:i- 


L 


memo  AMS 


nun 


(B.6) 


i-i  -  f, V 

o(-m) 


then  (B.3)  can  be  written  as 


K  ill  ,  i 

tL_  (_£)  p  (sin  e)  - — r  ("l)  [cos  m(t>  +  i  sin  mO) ] 

(2-8  ^ 

on 


(B.7) 


+  J,w  sLcos  o(l>  -  i  sin 

t(  -m)  ( 


Therefore  for  even  ro  (B.7)  becomes 


R  il-l 

^  P,,  (sin  f))' 

R  „  r  om 


L  (2-6  ) 

OM 


COS  m4)+i 


^  Im  ^  -m) 

I-  (“-'u)*  -* 


sin  (B.8) 

even  m 


and  for  odd  m 


R  (5+1 


^■g(  1 

(2-6  )M 

'  om  J 


cos  m$-i 


(2-6  )®  J 

om  -* 


sin  !»)> 


( ' Odd  m 

' 


(B.9) 


Then  for  even  in  let 


’  Cm 


even  ra 


(2-6  )*  ’ even  m 

'  om 


(B.IO) 


S,; 

Cm 


even  m 


(B.ll) 


even  m 
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and  for  odd  m  lot 


'£m 


odd  m 


(2-f>  Odd  ra 

on 


(B.12) 


(B.13) 


Substituting  (B.lO)  throuRh  (B.13)  into  (B.8)  and  (B.9)  and  summiiifi:  over 
the  required  positive  values  ol*  /  and  m  ^ives 


U(r,0,<t))  = 


0  >0=0  in=i 


Ul 


S  yj  sin  m(i> 

iJm 


(B.14) 


which  is  identical  to  (B.l). 

To  solve  for  J,  and  J,.  ^  in  terms  of  C,  and  S,  first  write 

im  f(-ni)  p.in  Pm 

from  (B.lO)  and  (B.ll)  for  even  m 


( 2  “  8  )  C  s  —  J  ,j  'I'  tJ  \ 

om  Zm  om  £(  -m) 


(B. 15) 


1(2  -  br  s„  =  j„  -  j 


on 


£m 


Im  ^(-m) 


(B.16) 


Adding  (B.15)  and  (B.16)  yields 

“im  "  -  ''om>^"'im  '  *  ’  ">■”> 

Equation  (B.17)  is,  however,  inaccurate  for  the  case  m  O  lor  which 
the  following  relation  holds 
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(B.18) 


To  incorporate  (B.18)  Into  (B.17),  write  2  as 


,  si  1 2  m  /  0 

(2  -  (2  -  )  =  j 

om  om  I , 

I  1  ra  =  0 


Then 


(2-8 


r?  i 

om^  ® 


ven  m 


Alter  subtracting  (U.16)  from  (B.15)  this  same  type  analysis 


oven  III 


' ‘  777711  I 

ora"^ 


even  m 


UsliiK  (n.l2)  and  (B.IB)  to  analyze  the  odd  m  case  yields 


om' 


odd  tn 


'  77771 

^2  -  'odd  m 


Rewriting  (B.20),  (0.21).  (B.22).  and  (B.23)  in  the  following 
pact  form 


(2  .  r,  )J  ‘  > 

om 
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(B.19) 

(B.20) 

yields  for 

.  (B.21) 

(B.22) 

(B.23) 

more  com- 

(B  24) 


It  Is  also  possible  to  re— express  as 


'  om 


(B.2G) 


«(-m) 


- — T  V"*” 


-iniX 


(2  -  ;.  ) 
'  om' 


(1!.27) 


From  (B.24)  and  (B.25)  the  following  equation  may  be  written 


-* 


m=-^ 


im  im 


-  Ha  S 

m=0 


(B.2«) 


The  expected  value  of  (B.28)  is  defined  as  the  niis  potential  coef  ("icient 
variation.  Therefore 


(It  .  2!)) 


This  gives  us  the  first  part  of  (3.15).  To  derive  the  Heller  model 
2 

counterpart  to  CT?  begin  by  specializing  (A. 11)  to  Heller's  three  laver 

I 

model 


U( 


"n-  “S.  -  ra  E  g  I 


(2<5-I'1)p  eos  V) 


('■+1 


(B. 30) 


Substituting  (A. 10)  into  (B.30) 


X  .  (B.31) 


Afcaiti,  Koinp:  back  to  App.  A,  Eq.  (A. 8) 

I) 

In  conjunction  with  (B.3i)  yields 


t. 

n 


II 


A'  ^R|-DA^*+l 


i6r  ^  J  '  %  f 


(B.33) 


But  recalling  the  identity 


I 


^TTb  3|  I 

i'  £  m*m 


(B.34) 


(B,33)  simplifies  to 


t  t(-p- 


7 


(B.35) 


X  U(l.„-D,,  O',  . 

Since  V  ^  and  m'  r  m  in  (B.35),  then  for  a  Riven  P,  and  tn  only 
one  term  In  the  r  ^  sequence  will  be  non'ero.  Therefore  (B.35)  may 
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be  written 


3  r 

-  h  eJ  ' 

:',  1  =  1  ■'  '  ' 


Y  ,  ,u;  i{  -i).,  ui.  . 

i.m  ’  i’  ’ 


Multiplying'  botli  sides  of  and  tlioii-  complex  coriiiinalcs  and  t 

the  expected  value  fiives 


■d-  r 


1  1 

3 

V  t 

j  2 

lu 

^  T 

y  \  1^,  > 

(tbir 

i=l 

i'=l 

•  I  1 


Im 


...i 


X  U(|{^-D.,  s  ■:)  cii-  U3.|. 


Because  the  potential  variations  fi-oni  different  slu'lls  arc  uiicorrcl 
then 

I 


E  U{R-D^,  3')  U{i{  -1)^,  U  '  l\  '  U  'd-  '  I 


I 


fi>r  i  /:  i' 
lor  i  i' 


where  is  any  function  of  '  and  '  .  Apnlyinu  (B.B8)  to  (B 


- )  CT- 

U  /  .1  . 


t  ^ 

i>,ir  id  \  /  111  -1). 


Y  •  (  ,  :  'Y  V--,  :  ;cj 

m  ’  m 


However,  the  inte^i’nl  in  (IJ.Hf))  <'qiia1s  -I".  sim  pit  fxinu 

usiu'^  tliis  fact  r:>ves  the  desirtMl  result 


2 


(n.y,  3  /b-i)x.h-. 


(It.  :t(  1 


aki  n;;' 


(B,;t7) 


a  I  ed  . 


(B.  ;tK) 


.  B7) 


(B  .  :th) 


(15  Bh) 


(15.  10) 
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Appendix  C 


LOW  FREQUENCY  FILTER  CONSTANT  AND  RMS  POSITION  ERROR 
ESTIMATE  ERROR  DERIVATIONS 


The  measurement  combinations  used  in  tills  research  are  ^iven  below  by 
number  preceded  by  the  assumed  direction  of  estimation.  Note  that  these 
derivations  assume  the  use  of  the  Heller  model. 

CT-1)  If  a  measurement  of 


z  = 


(C.l) 


Is  taken  then 


and 


H  =  [0  1  00  0] 


(C.2) 


R 


(C.3) 


Equation  (4.70)  gives  the  root 


a 


Using  (4.76)  to  determine  <X 


then 


(I)(_4k>-..i^2k^  =  0 
111 

s 


(C.4) 


(C.5) 


(C  .(  ) 


To  obtain 


(fr~  )  the  following  integral  (from  4.49)  must  be  evaluated 
^y  P. 


t  T  T 

K  ^ 

y 

- —  fis 

(-) 


1 

2n.1 


B 

2k  sfy 

— -  ds  . 

W^Cs  4  fr) 


(C.7) 


Usinp  the  Cauchy  residue  theorem  this  equals 


1 

2ni 


2kgSa 


0)  (s  4  ■^) 

s 


ds 


2k®0<p. 


2 

0) 

s 


(C.8) 


or 


(C.9) 


To  Include  the  EVM  in  this  and  any  of  the  subsequent  total  estimate  error 
variance  formulas  the  only  modification  necessary  Is  to  change  k  to  k'. 


CT»«2;  In  the  case  where 


then 


and 


“o  1  0  0  0 
0  0  0  1  0 


R 


(C.IO) 


(C.ll) 


(C.12) 
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The  low  frequency  pole,  obtained  from  (4,70),  Is 


f-  = 


A  A  A 

U4  8k  k,.  -  :i6k. 
3  4 


- ^,1 


(C.13) 


For  multiple  measurements  the  vector  O  Is  obtained  from  both  (4.76) 
and  (4.77).  Equation  (4.76)  fjives 


- 1 

1 

to 

o  o 

_ 1 

■«r 

1 

2 

U) 

■-r 

0 

L.  .J 

L  2j 

s 

(f'.14) 


then 


O'. 


To  determine  use  (4.77) 


wA  ^  2 

•‘a  “"s 


ri  -  2k^s2  6kjs 


A  A 

-6k^s  24k5  +  r^^ 


a 


k. 


-  0. 


s=— 


(C.15) 


(<M6) 


The  second  row  of  this  equation  Rives 


0:„ 


6k^-flV 
4  X 


^  *'2^ 


(C.17) 


The  Inst  tntoRral  In  (4.40)  <'an  now  be  «'xpresso<l  ns 


. 
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1 

2n.i 


t  T  T 

9  (H)ll  (h) 

K  r*  * 

y 

-  ds 

(-) 


1 

2n.j 


w  (s  +  -fi) 
s 


ds 


The  evaluation  of  (C.18)  gives 


+  6k 

2 

U) 

s 


(C,19) 


IT-1;  If  the  measurement  Is 


=  r  +  e 

XX 


(C.20) 


Then 


H  =  [1  0  0  0  O] 


(C.21) 


and 


(C.22) 


Tho  characteristic  equation  from  (4.70)  Is 

kjS  4  r^  =  0  . 


(C.23) 


Solving  for  the  two  left-half  plane  poles  results  In 


<jR+) 


(C.24) 


where  L  is  a  low  frequency  factor  coefficient  defined  as 


1. 
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(C.25) 


■  'w^wnifvwr...^ «... 


To  determine  d  and  e  write  (4.79)  as 


ki  s 


2  2 
s  +  T’S 

Li 


(ds  +  e) 


*larfce  s 


2  '"l 

U) 

s 


lar^e  s 


(C.26) 


Using  the  approximation 


2  2 
s  +  — s 

Li 


s  Ls 


(C,27) 


(C.26)  may  be  written 


small 


(C.28) 


Zeroing  the  first  and  second  degree  s  terms  gives 


d  = 


w 


(C.29) 


2V 


e  = 


L<i' 


(C.30) 


The  low  frequency  rms  variance  is  obtained  by  integrating 


9^  (s)hV,(s) 


ds  =  + 


(-) 


J-  ^  ^ 


n  3 

k,s  (ds  4  e) 


(0^(  >4  ) 
S 


ds.  (C,31) 


Summing  the  residues  at  both  lel't-half  plane  poles  in  (C ,  31 )  gives 


Zk^e 

I  (0^1.^ 


(r.32) 
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11=2-.  a  incasiiremont  of 

1 


then 


anfl 


r  -r  )4-  fr  ^IrVe, 

2  XX  yy  '  xx  2  zz  ’ 


H  r  I  O  0  0  1 


H  . 


(C.33) 


(C.34) 


(C.35) 


TIh'  chafaftorlstlc  equation  from  (4.70)  Is 


A  4  A  2  /  A 

‘‘i  ”  ■*  ■'  \^^ry  '‘l 


0 


(C.36) 


The  left-hair  plane  poles  of  (C,34)  produce  a  low  frequency  factor 
expressible  ns 


(f.H) 


2.a^6i. 


S  +  U)., 


(C.37) 


where 


k.  r 


A  2 


5  1±  (A 

I  1  1 


,i,(i 


1/2 


\  a)  ’  5 

L  “l 


1/4 


(C.38) 


(C .  39) 


(C.40) 


The  convergence  condition  for  this  situation  is 
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SummlnK  the  residues  at  both  of  the  left-lialf  plane  poles  in  (C.44) 
ultimately  plves 
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(c  VI?) 


ts. );  ■  [ 


4  k  Oj 


4  (8j-ft^)]  -2k'^l58j  ^ 


where 


co«(i 


(C.48) 


IT- 3:  h'or  a  moasiiretnent  of 


1 - 

(- 

1 

(- 

r  4  -^r 

^  XX  yy 

+ 

XX  2  zz 

-?(r  -  r  ) 

-i?  4  r 

2  7.Z  yy  _ 

2  XX  zz 

4  e 


(C.49) 


t  lie  n 


I  0  0  0  I 
i  o  0  0  1 


(C,50) 


It  ^ 


*1  " 


^oJ 


(C.51) 


The  eharaetoristlc  equation,  aifaln  from  (4.70),  reduces  to 
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where  It  has  lieon  assumed  that 
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If  Tj  Is  small  then  the  s  term  in  (C.52)  is  ncKllKible  and  the 
equation  becomes 
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Solving  for  the  left— half  plane  roots  the  low  frequency  factor  is 
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where 
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Here  the  convergence  condition  (4.79)  is 
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SubstitutiiiK  (C.27)  Into  (C.58)  yiclfls 
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UsiiiK  only  the  first  equation  In  (C.59)  since  the  two  are  linearly  de¬ 
pendent  (jives 
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7eroin(;  the  first  and  second  de(jree  terms  in  (C.60)  gives 
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b) 


This  results  in  four  unknowns  but  only  two  equations.  The  remaining  two 
equations  come  from  (-1.80)  expressed  ns 
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For  small  s  and  small  (C.63)  rediucs  to 
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Using  only  the  first  of  the  equations  in  (('.(14)  because  of  their  linear 
dependence  gives 
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This  immediately  produces  the  conditions 
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Solving  (C.61),  (C,62),  (C.66)  and  (C,67)  simultaneously  produces 
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The  analytic  low  frequency  variance  integral  reduces  to 
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Solving  (('.72)  again  with  the  aid  of  Cauchy’s  theorem  the  analytic 
variance  equals 


\  x/jf  w  \l,  L  ' 


(C.77) 


-122- 


HEKEHKNCKS 


1.  Hildebrant,  R.R.,  K.R..  aiul  MinRlfii ,  S..T..  "  I'lie  Ellccls 

of  Gravitational  lincerlaintios  on  t  lie  Ertors  ol  I  ortial  Navi^cation 
Systems",  J.  of  the  Inslit'ite  of  Navigation.  Vol  .  21,  No.  4, 

Winter  1974-1975. 

2.  Bell  Aerospace  Corp.,  "Gravity  G radiometer  Instriiment  ieasihility 
Demonstration  and  Studies  and  Gravity  Gradiometer  System  Definition 
Studies  ProKram",  Bell  Aerospace  Corp.  ,  Niattara  l  alls,  NY,  Vol.  1 
Tech.  Kept,  No.  D6441-953001 ,  Oct.  1976. 

3.  "Advances  in  Dynamic  Gravimetry",  Proc.  of  the  Symposium  on  Dynamic 
Gravimetry.  Mar.  16-17,  1970,  Kort  Worth,  Texas. 

4.  Pelka,  E.J.,  "On-Line  Parameter  Estimation  and  Control  for  a  Rotat¬ 
ing  Gravity  Gradiometer",  Ph.D.  Dissertation,  Stanfortl  University, 
Dept.  Aeronautics  and  Astronautics,  Guidance  &  Control  I.ab.  . 

Aug.  1977. 

5.  DeBra,  D.B.,  Schaechter,  D, ,  and  K'lrosaki ,  M..  "Study  to  Develop 
Gradiometer  Techniques",  Air  Force  Geophysics  Laboratory  Technical 
Report,  AFGL-TR-78-0046,  Dec.  1977  (AFGL/LWG,  Hanscom  AFB,  MA  01731), 
December  1977. 

6.  DeBra,  D.B.,  Breakwell,  J.V.  ,  and  Kurosaki ,  M. ,  "Study  to  Develop 
Gradiometer  Techniques,"  Air  Force  Geophysics  Laboratory  Technical 
Report,  AFGL-TR- 79-0063,  January  1979. 

7.  Heller,  W.G,,  "Free-Inertial  and  Damped- Inert! al  Navigation  Mechan¬ 

ization  and  Error  Equations",  The  Analytic  Sclcniccs  ('orp.  ,  Reading, 
Mass.,  TR-312-1-1,  18  Apr.  1975. 

8.  Rice,  D.A,,  "A  Geoidal  Section  in  the  United  Slates",  Bulletin 

Geod. ,  Vol.  65,  1962. 

9.  Kasper,  J.F.,  Tr.,  "A  .Secoiul  Order  Markov  Gravity  Anomaly  Model", 

J.  of  Geophysical  Researcli,  Vol.  76,  No.  32,  Nov.  1971. 

10.  Heller,  W.G.  ,  "Gradiometer— Aided  Inert  ial  Navigation",  The  Analytic- 

Sciences  Corp.,  Reading,  Mass.,  TR-312-5,  7  Apr.  1975. 

11,  Heller,  W.G.,  "A  New  Self-Consistent  Gtalistical  (iravity  Field 
Model",  Proc,  of  the  American  Geopnysicral  Union.  Fall  Annual 
Meeting,  San  Francisco,  CA.,  Dec.  1976. 

12,  Breakwell,  J.V,,  "Satellite  Determination  of  Short  Wavelength 
Gravity  Variations",  J.  of  Astronaut ica 1  Sciences,  Vol.  37,  No,  4, 
pp,  329-344,  Oct-Dc'c.,  1979. 

13.  Breakwell,  J.V,,  Kurosaki,  M.,  "Real-Time  Gravity  Filtering  I'rom 
On— Board  Gradiometers ,  "  AIAA  GttiC  Coni'.,  Boulder,  Colorado, 

Aug.  1979, 


-12.3- 


i 

j 

I 


14,  Lorch,  I'..!.,  Klo.sho,  S.M.,  Laubscher,  U,E.,  and  Wagner,  C.A,, 
"Gi'avlty  Model  I  nproveinciit  Using  GEOS— 3  (GEM  9  &  10)",  Goddard 
Siiacc  Kliglit  I'entoi’,  Greenbelt,  Md .  ,  Scpi,  1977. 

15.  Ttagcser,  M.B.,  "A  Status  Report  on  the  Floated  Gravity  Gradiometer 
lU'pL.  No.  (’  3868,  The  C.  Stark  Draper  Iiaboratory,  Cambridge,  Mass. 
Nov.  30,  1972. 

16.  Bryson,  A.E.,  and  llo,  Y.C.,  Applied  Optimal  Control.  John  Wiley  & 
Sons,  Itie.,  N('W  York,  NY,  1975. 

17,  Iloiskancn,  W.A.,  and  Moritz,  U.  ,  Pliysical  Geodesy.  W.H,  Freeman, 
San  Francisco,  CA.,  1967. 


